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CHAPTER  IV 


D I R E C T L Y- T R A N S M I T T E D  SROUND  PHENOMENA 
ELASTIC  ANALYSIS 

A.  INTRODUCTION 

This  ebsptST  consists  of  the  report  by  N.  M.  Nevmsrk  sod  Associates  of  the 
results  of  a  suboontraot  study  under  Cimtraot  IlA-49~146-XZ-073«  Only  minor 
editorial  changes  hare  been  made  to  fit  this  stu^  into  the  format  of  the  overall 
final  report  under  this  contract*  The  report  vas  prepared  by  N.  M.  Neemark  end 
A.  Ang  vith  the  assistance  of  J.  F.  Martha  (^o  contribtxted  the  section  on  the 
equation-of-state  effects),  A.  R.  Robinson  (idio  contributed  the  analytical  solu¬ 
tion  of  the  spherical  problem),  and  S«  Sutcliffe  (i^o  contributed  the  section  on 
the  spherical  solid  vith  bilinear  behavior),  Coioputer  progranming  for  these 
studies  vas  done  by  J.  V.  Melin,  G.  N.  Harper,  and  J,  Rainer, 

1,  Statement  of  Rroblem  and  Objectives 

The  prediction  of  close  range  ground  motions  and  pressures  resulting  from 
nuclear  explosions  involves  a  vide  range  of  material  properties,  ranging  from 
that  of  a  purely  liquid  state  to  that  ^ich  is  essentially  of  the  original  solid 
state.  The  transition  betveen  these  tvo  states  of  the  material  is  gradual  and 
probably  no  definable  interfaoe  exists  betveen  any  tvo  intermediate  states. 

For  purposes  of  studying  the  effects  of  ground  shocks,  the  total  earth  motions 
generated  by  a  nuclear  blast  are  iisually  divided  into  those  that  axe  caused  by 
stress  vaves  vhich  are  directly  transmitted  throu^  the  earth  material  from  the 
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•ntrgy  r«l«M«d  dir«eily  into  th*  •arth  by  a  blast,  and  thosa  that  ara  induead 
by  atrass  vavas  idiloh  ara  initiatad  from  tha  aarth  swfaoa  by  tha  axpanding  air 
shooka. 

Baaad  on  a  hydrodsmamio  modal,  Eroda  and  Bjork*  hava  Buooaasfully  performad  numari 
oal  oaloulationa  for  a  simulatad  two  magaton  burst  in  tuff,  a  soft  rock  material^ 
Erode  and  Bjork  pointed  out  that  'tiiair  calculations  for  tha  low  prassura  rangas, 
lass  than  8  kilobars,  ara  not  strictly  -valid  sines  in  -thssa  prassura  rangas,  -tha 
affaets  of  -fcha  plastic  and  elastic  properties  of  -the  solid  material  become  signifi 
cant. 

Using  -the  results  of  Erode  and  Bjork  in  -the  transition  pressure  range,  (assumed 
to  be  at  a  radial  distance  of  660  feet  from  the  center  of  a  2  megaton  burst),  as 
a  s-tarting  condition,  directly  induced  ground  motions  for  a  half-4pace  solid  were 
determined  numerically  on  the  assumption  that  -the  material  beyond  the  region  of 
the  crater  is  homogeneous,  isotropic  and  linearly  elastic.  Such  an  assumption, 
of  course,  is  not  realistic  since  a  plastic  region  would  precede  -the  elastic 
region)  howe-ver,  the  use  of  this  assumption  was  dictated  by  a-vailable  means  for 
ob-taining  approximate  answers.  Results  of  an  approximate  study  of  -the  effect  of 
-the  equation  of  s-tate  indica-ted  that  no  significant  error  in  -the  applied  pres¬ 
sures  is  in-vol-ved  by  neglecting  a  plastic  region  between  the  liquid  and  the 
elastic  regions.  Aside  from  -the  half-space  solid,  -the  simpler  problem  of  an 
infinite  space  subjected  to  an  explosi-vw  pressure  applied  in  a  spherical  ca-vity 
was  sol-ved. 

In  this  study,  no  consideration  of  the  effects  of  -the  expanding  air  blast  pres¬ 
sures  on  the  surface  was  inclvided  althou^  this  can  be  considered  without  much 


*Brode,  H.  L.  and  Bjork,  R.  L.,  "Cratering  Fr<Hn  A  Megaton  Surface  Explosion", 
RM-2600,  The  RAND  Corp.,  June  30,  1960 
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difficulty.  The  tolution  of  on*  of  tho  bolfHipao*  probloma  izrrolTOS  a  oituatioa 
idiioh  is  osBontially  similar  to  that  of  tho  actual  air  blast. 

Tho  spocific  inrobloms  of  vhioh  solutions  aro  prosontod  horo  aro  outlinod  bolotrt 

Problom  1 — A  full-spaco  vith  a  spherical  hole  stibjectod  tc  a  uniformly  distri¬ 
buted  pressure  (see  Figure  9te). 

Rroblem  2— A  half-space  vith  a  stress-free  semi-spherical  oavityi  loads  are 
applied  on  the  horisontal  surface  ^ioh  are  exactly  equal  to  tho  negative  of  the 
tangential  stresses  from  Problem  1,  referred  to  as  "correction"  loading  (as  shown 
in  Figure  98b).  Resvilts  from  this  problem  aro  referred  to  in  this  report  as 
"correction"  stresses,  velocities,  or  accelerations. 

Problem  3--A  half-space  vith  a  semi-spherical  cavity  subjected  to  a  \iniformly 
distributed  blast  pressure  applied  directly  on  the  boundary  of  the  cavity  (see 
Figure  98c). 

Problem  4— The  same  as  Problem  3,  except  that  the  applied  pressure  is  uniformly 
distributed  only  over  the  lover  half  of  ibe  semi-spherical  surface  of  the  cavity; 
the  pressures  on  the  upper  half  of  the  cavity  are  distributed  as  p^  sin  20  (as 
shown  in  Figure  98d). 

2.  Methods  of  Analysis  and  Results 

Problem  1  vas  formulated  analytically  and  caloulations  for  the  radial  and  tangential 
stresses  vere  performed.  Additional  results  vere  also  obtained  by  an  alternate^  analytical 
procedure  for  the  infinite  space  subjected  to  an  eiqponentially  decaying  pressure¬ 
time  curve  (see  Appendix  1).  Solutions  for  stresses  in  IToblem  1  obtained  vith 
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th*  uaa  of  a  siiharloally  ■ymnatrlo  diserote  modal  vara  in  vary  good  agraamant  vith 
tha  oorrasponding  analytical  aolutiona. 

Rasults  for  the  half-apaoa  aolids,  Problama  2,  3  and  4,  vara  obtained  axoluaivaly 
by  tha  technique  of  digital  simulation.  Although  no  analytical  solutions  are  avail¬ 
able  for  direct  verification  of  these  results,  there  are  a  number  of  limiting  con¬ 
ditions  ^ich  all  the  solutions  satisfied,  thus  indicating  the  correctness  of  the 
solutions.  Furthermore,  the  mathematical  consistencies  of  the  models  provide  a 
basis  for  ascertaining  the  correctness  of  the  solutions  if  the  numerical  results  cure 
reasonable  from  an  intuitive  standpoint. 

All  calculations  were  performed  on  large  hig^-speed  digital  coiiq>uters.  In  addition, 
>dienever  the  discrete  models  cure  involved,  the  problems  were  simulated  directly  on 
the  computers. 

All  results  are  presented  graphically.  For  the  full-space  solid,  time  functions  of 
the  radial  and  tangential  stresses,  and  of  the  particle  displacements,  velocities, 
and  accelerations  at  points  of  varying  radial  distances  from  the  burst  are  given. 
Analytically  determined  solutions  cure  given  only  for  the  rctdial  and  tangential  stresses. 
For  the  half-space  solids,  the  radial  cuid  tcuigential  stresses,  and  the  radial  pcurticle 
motions  are  plotted  against  time.  These  vere  presented  for  points  along  lines  vith 
different  angular  positions  from  the  horizontal.  The  solutions  for  the  half-spaces 
include  other  stresses  (shear  and  circumferential)  cmd  the  tangential  accelerations, 
velocities,  and  displacements;  hovever,  these  are  not  presented  here. 

3.  Stcurting  Conditions 

The  original  data  from  the  calculations  of  Brode  cuid  Bjork,  previously  mentioned, 
vere  presented  in  the  form  of  contours  connecting  points  of  equal  pressure. 


Ttlooiiy  or  doniity  at  apaoifitd  times*  From  these  oontours,  pressure-time 
relationships  were  reproduced  and  presented  in  a  form  that  is  more  meaningful 
for  further  study  and  analysis.  These  pressure-time  relationships  at  radial 
distances  of  330,  495,  and  660  feet  from  the  point  of  bvirst,  are  presented  in 
Figures  99  throu^  101.  Corresponding  relationships  for  the  particle  velocities 
were  also  reproduced  but  are  not  presented  here. 

The  problem  considered  in  this  study  involve  a  full  space  with  a  spherical  hole, 
or  a  half-space  with  a  semi-spherical  cavity  on  the  surface.  The  radius  of  the 
cavity  in  all  cases  is  660  feet  from  the  center  of  a  tvo-megaton  burst. 

The  jtressure-time  relationship  of  the  applied  pressure  is  shown  in  Figure  102, 
%diich  represent  the  pressure-time  curve  from  the  calculations  of  Brode  and  Bjork 
at  660  feet  directly  below  the  burst.  The  decay  i>ortion  of  this  curve  is  extra¬ 
polated  by  inference  from  available  measured  data.  In  the  solutions  by  digital 
situilation,  this  pressure  is  actvially  applied  at  a  radial  distance  of  645  feet 
from  the  point  of  b\irst,  and  not  at  660  feet,  for  convenience  in  treatment  of  the 
lumpedHsass  model.  Ihere  is  therefore  a  sli^t  discrepancy  between  the  analytdcal 
and  the  lumped-mass  solutions;  this  is  not  of  significance  in  the  results. 

The  velocities  at  the  surface  of  the  cavity,  ^ich  were  also  calculated  by  Brode 
and  Bjork,  were  not  considered  in  this  study, 

B.  EFFECTS  OF  THE  EQUATION  OF  STATE 

In  the  present  study  of  the  stress  field  outside  the  crater  region,  the  pressure¬ 
time  relations  at  some  radius  as  obtained  from  the  hydrodynamic  model  of  Brode  and 
Bjork  are  used  as  the  boundary  stresses  applied  to  an  elastic  solid.  Theor¬ 
etically,  the  pressures  derived  from  the  hydrodynamic  model  results  should  be 

adjusted  in  order  to  provide  a  consistent  boundary  condition.  The  need  for  this 
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adjuBtownt  rtiultt  from  two  f motors:  nsmsly,  (l)  tho  wsto  propagation  valoolty 
in  tho  tlastio  solid  vill  bo  loss  than  that  in  tho  fluid  modal ,  and  (2)  tho  par- 
tiolo  volooity  at  tho  boundaiy  botwson  tho  fliiid  and  solid  rogions  should  bo  tho 
samo  in  both  media*  That  thoso  tvo  conditions  aro  not  autcnnatioally  satisfied 
is  tho  result  of  omitting  the  plastic  or  transition  region  idiioh  in  tho  real 
situation,  separates  tho  hydrodynamic  and  elastic  regions*  It  is  apparent  that 
any  procedure  adopted  to  adjvist  these  pressures  can  only  be  an  approximation* 

In  the  subsequent  discussion  it  is  shovn  that  the  necessary  adjustment  appears 
to  be  small  and  may  be  neglected  vithout  serious  error* 


If  the  mathematical  model  is  assumed  to  consist  of  compressible  fluid  and  elastic 
solid  regions,  a  vave  propagating  through  the  fluid  could  be  partially  reflected 
at  the  boundary  separating  the  tvo  regions*  The  relations  idiich  must  be  satis¬ 
fied  by  the  incident  fluid  pressure  and  particle  Telocity,  p  and  and  the 
transmitted  stress  and  particle  Telocity,  a  and  d,  are  as  follows: 


p  +  Ap  =  a 

(la) 

II 

(lb) 

idiere  Ap  and  Ad^  are  changes  in  the  fluid  quantities  resulting  txom  the  reflec¬ 
tion*  The  relationships  between  stress  and  particle  Telocity  in  both  the 
incident  fluid  and  solid  waTes  are  time  dependent;  howeTer,  at  the  front  of  an 


ideal  wave  the  following  relations  exist 

d  = 


a 

'  8  8 


~  P-C 


(2a) 

(2b) 


f~f 


*  Force  equal  change  in  momentum  per  unit  time. 
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iih«r«  p  ud  e  ar«  dtaaity  tad  imfv*  TvloelVi  mptetlTvly,  and  iht  nibteripts 
rafar  to  tha  aolid,  £,  and  fluid,  f.  Althoui^  tha  fluid  praaaura  pulaa  ahnpaa 
oaloiilatad  by  Broda  and  Bjork  ara  not  truly  idaal,  Equation  (2b)  ia  approiciiBataly 
oorraet  for  tha  initial  part  of  tha  pulaa.  Tha  application  of  Equation  (2b)  im- 
pliaa  a  Talua  of  tha  paraiaatar  p/ft^  ^ioh  is  not  time  dapandant  for  a  given 
radius.  Tha  folloving  tabla  suomarisaa  tha  valuaa  of  this  paraaatar  for  tha 
prassura-tima  and  valoei-fy-tima  pulsas  in  tha  fluid  at  a  ranga  of  200  matars: 


Tima 

m  sac 

P 

kilobars 

(kilobar-aao  )/kiloiDatar 

58 

3.4 

45.4 

65 

5.0 

48.6 

74 

4.5 

45.0 

82 

3.0 

32.0 

85 

2.4 

31.5 

Thaaa  data  indioata  that  tha  assumption  of  a  constant  value  for  p/d^  is  raason- 
abla,  at  laaat  for  tha  principal  part  of  tha  praasura  pulsa. 


If  it  is  assimad  that  a  ralation  similar  to  Equation  (2b)  is  valid  for  tha 
raflaotad  praasura  and  valoci'fy,  i.a., 


^fr 


(3) 


tha  ralation  batwaan  tha  incidant  and  transmittad  strassas  is  as  follows: 


whara 


a 


-  _Se_ 

■  1  + 


(4) 
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In  Equation  (4)  a  repirosents  the  transmitted  stress  that  is  dsTsloped  vfaen  a  plane 
vave  impinges  on  the  boundary  betveen  tvo  different  media.  As  such,  it  would  be 
correct  at  least  at  the  shock  front. 


Since  Equation  (4)  represents  an  acoustic  reflection,  it  Is  desirable  to  calculate 
the  reflection  factor  for  plane,  finite  soqplitude  waves,  in  the  fluid  reflecting  from 
a  rigid  vail.  The  equation  of  state  used  by  Erode  and  Bjork  for  tviff  is 

p  =  (0.425E  +  0.113n^/^)  +  5.30n  10^° 

10-^^ 

0.707n  (■^)^:  10^° 

^ - 10“  ^5) 

10^  +  (E/10^°) 

^diere  p  =  pressure  in  dynes/cm^ 

E  =  specific  internal  epergy  in  ergs/gm 

p 

n  =  !  p  =  initial  density  =  1.7  gm/cc 

fo  ° 

For  pressures  less  than  10  kilobars  and  n  less  than  about  1.3,  Equation  (3)  may  be 


approximated  as 


5.3  , 


Combining  Equation  (6)  with  the  Hugoniot  energy  equation. 


J  (P  -  P„)  3 


results  in  the  following  pressure-density  relation: 


p  =  8.23  X  10 


P  2  P 


dynes 

2 

cm 


for  the  ambient  pressure,  p^,  equal  to  zero.  As  a  check  on  Equation  (8),  we  may 
calculate  the  wave  velocity,  c^,  for  waves  of  small  aiiq>litude; 


7.200  ft/see 


c 


0 


for  yory  small  values  of  —  .  The  value  of  usually  quoted  for  tuff  ranges 
from  about  6,000  -  8,000  ft/sec. 


When  a  finite  amplitude  vave  undergoes  normal  reflection  from  a  rigid  boundary, 
the  ratio  of  the  excess  pressures  behind  the  reflected  and  incident  vaves  is  given 
by  Cole*  as  follovs  ^ 


f 


1  + 


^  _ i 


(9) 


vhere 


is  the  density  behind  the  reflected  shock  front. 


The  ratio  of  reflected  to  incident  pressures,  p^/p  may  be  calculated  using 

Equations  (8)  and  (9).  For  an  incident  pressure,  p  =  9  kilobars,  p  /p  =  2.2. 

r 

For  the  same  incident  pressure  in  vater,  p^/p  =  2.9.  Thus  the  reflection  factors 
for  finite  amplitude  waves  in  the  range  of  interest  do  not  appear  to  greatly  exceed 
the  acoustic  reflection  factor,  idiich  is  2.0. 


Approximate  values  for  the  wave  velocity  in  the  fluid  can  be  derived  in  several 
ways,  each  of  which  yields  about  the  same  result.  Using  ideal  shock  relations 
and  the  calculated  maximum  values  of  p  and  u^  from  Brode  and  Bjork  the  following 
results  are  obtained: 


Range  f^ 

Meters 

150 
200 
250 

*Cole,  R.H.,  "Underwater  Explosions' 

New  Jersey,  1948. 


ft/sec 
8,650 
8,100 
8,030 

Princeton  Ikiiversity  Press,  Princeton 
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Thus  at  a  range  of  200  meters  the  appropriate  value  of  =  8,100/7»200  =  1.12 
and  from  Equation  (4)  o  =  0.945  p.  This  reduction  in  the  fluid  pressure  is 
sufficiently  small  that  it  may  be  neglected.  Thus  the  applied  stress  for  the 
elastic  solid  may  be  taken  as  the  fluid  pressure  at  the  radius  of  the  boundary. 

The  applicability  of  the  results  of  Erode  and  Bjork  as  boundary  conditions  for 
subsequent  oalotilations  in  media  other  then  tuff  depends  upon  several  factors 
including  the  effect  of  differences  in  the  equation  of  state  between  any  other 
material  and  that  assumed  for  tuff.  Of  partiouler  io^rtanoe  is  the  effect  of 
the  equation  of  state  on  the  pressure^ime  variation  from  point  to  point. 
Unfortunately,  no  calculations  are  available  for  any  other  rook  or  soil  material.  A 
qualitative  estimate  of  the  effect  of  material  properties  can  be  made  by  oonqparing 
the  peak  overpressure-distance  relation  for  water  and  a  gas  with  the  results  for 
tuff. 

The  variation  of  peak  pressure  with  vertical  distance  below  ground  zero  is  given 
for  tuff  in  Figure  103  along  with  similar  results  for  water  and  an  ideal  gas  with 
^=3.  The  RAND  Corporation  calculations  for  tuff  were  terminated  at  a  peak 
pressure  of  about  5  kilobars.  At  these  lower  pressures,  the  peak  pressure  varies 
approximately  as  the  inverse  three-^lves  power  of  radial  distance,  ndiile  at 
hi^er  pressures,  not  shown  in  the  figure,  the  peak  pressure  variation  is 
approximately  like  the  inverse  cube  of  the  distance.  The  dashed  lines  in  the 
figure  are  extrapolations  of  available  data. 

The  ideal  gas  solution  was  also  obtained  by  Brode  and  Bjork  using  the  hydrodynamic 
model.  Here  the  peak  presswe  varies  as  approximately  the  inverse  ovd>e  of  the 
distance  for  all  the  ranges  calculated.  A  comparison  of  the  tuff  and  ideal  gas 
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oaloul»tion«  iodimt*!  that  for  •  giYon  roogo  the  peek  preeeuree  between  the  two 
meter i els  differ  by  less  then  e  feotor  of  ebout  three. 

The  reletlonshlp  for  weter  shown  in  Figure  103  wes  deri-ved  from  the  peek  weter 
OTerpress\]re*^i8tenoe  dete  for  e  deep  underwater  e3q;>lo8ion*.  In  deriving  the 
reletion  for  water  it  was  assumed  that  the  weapon  yield  was  1200  KT  or  twice 
the  total  energy  assigned  to  the  rock  half-spaoe  at  time  zero  by  Brode  and  Bjork. 

The  peak  water  pressure  -varies  approximately  as  the  inverse  of  the  distance  for 
pressures  of  about  0.1  kilobar.  As  in  the  previous  oasSf  less  than  a  factor  of 
three  separates  the  curves  for  -the  tuff  and  water  for  peak  pressures  ranging  from 
about  0.1  to  10  kilobars. 

On  -the  basis  of  these  comparisons  it  seems  reasonsble  to  es^ot  that  differences 
in  the  equation  of  state  for  -various  rocks  and  soils  would  not  cause  differences  in 
-the  peak  pressures  found  using  -the  hydrodynamic  model  of  more  -than  an  order  of 
magnitude  and  perhaps  much  less.  This  observation,  along  with  the  preceding 
discussion  of  the  boundary  conditions  for  the  elastic  model  suggest  an  approximate 
me-thod  of  extending  the  RAND  Corporation  calo\ilations  to  materials  other  than 
tuff.  Since  the  peak  pressure^istance  data  do  not  appear  to  be  too  sensiti-ve  to 
significant  variations  in  the  equation  of  state,  it  may  be  ass\>med  that  the  tuff 
res\ilts  of  Brode  and  Bjork  are  applicable  to  other  types  of  rook  and  perhaps  soil 
materials.  Eqiiation  (4)  may  be  iised  to  adjxist  the  pressure  input  for  the  elastic 
body  to  cooqtensate  for  the  differences  in  density  and  wave  velocity  between  tuff 
and  a  given  material.  It  should  be  enqhasized  -that  -this  procedure  is  recomnended  here 
because  no  close-in  calculations  are  available  for  materials  other  -than  tuff.  When 
such  calculations  ha-ve  been  made,  another  procedure  should  be  de-vised. 

^lasstone,  Samuel,  The  Effects  of  Nuclear  Veepona.  IT.  S.  Atomic  Energy 
Comnission,  April,  1962 
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In  all  of  the  probleoB,  the  material  of  the  solid  is  assumed  to  be  linearly 
elastio*  Values  of  the  elastic  oonstants  of  the  solid  are  assumed  to  be  ai^oxi- 
mately  those  of  materials  similar  to  tuff,  a  material  in  idxioh  the  oaloulations 
of  Brode  and  Bjork  apply<  These  Talues  vere  dsteimined  or  inferred  from  aTail- 
able  field  data,  and  are  as  follows: 

Lame's  constants:  %  =  3.9  x  10^  psi 

=  3.5  X  10^  psi 

Dilatational  wave  velocity,  c  =  6,000  fps. 

These  values  correspond  to  a  hnmg's  modulus  of  8.8  x  10^  psi  and  Poisson's 
ratio  of  V  =  0.265. 


C.  THEORETICAL  SOLOTION  FOR  INFINITE  SPACE 

Consider  an  infinite  homogeneous,  isotropic  elastic  space  vith  a  spherical  cavity 
of  radius  a.  Let  the  cavity  be  loaded  by  a  normal  stress  (0,  0,  t),  idiere 
0  and  0  are  angles  locating  a  point  .n  the  surface  of  the  sphere.  The  stresses 
will  lutually  be  compressive;  $  however,  the  ordinary  sign  convention  of  the  theory 
of  elasticity,  tensile  stresses  positive,  will  be  used  for  convenience.  In  idiat 
follows,  no  distinction  is  made  between  the  initial  and  final  configurations  of 
the  medium,  so  that  the  ordinary  theory  of  elasticity  applies. 


For  a  radially  symmetric  load,  it  is  apparent  that  the  only  displacement  is  the 
radial  displacement  u.  There  are  only  three  non-sero  stress  components,  a^, 
Ogg  and  C00.  For  small  displacements  the  acceleration  in  the  radial  direction 
is  a  =  ii,  the  strain  £  =  4^,  =  ~  .  The  equation  of  motion  in 

the  radial  direction  is 


1.x (A 

_2  Ar  '  rr'  r  ”  "r 


r2  ^ 


or 


2a^r  -  <^ee  -  ... 

+ - = - “  =(>\ 


dr  r 

Expressing  the  stresses  in  terms  of  the  strains,  ve  have 

»,,  =  »0i,  =  li  +  ^  (^T  *  r>  *  5^ 

The  eqmtion  of  motion,  Equation  (10),  becomes 


(10) 


{?<  +  7y)  +  2(  ^  +  ^)  i  -2(^  +^)^=pu 


or 


2  6  u  2u  /  \  ..  1  .. 

77*717-7- <TT^'”=^« 


(11) 


Nov  ve  introduce  the  displacement  potential  0  such  that  u  =  The  left  side  of 

dr 

Equation  (11)  becomes 

ar^  -  ar2'r2lt=  ar^^^2^r 

The  ri^t  side  is  just 

_i  X  (-i^) 

so  that  the  equation  of  motion  is  satisfied  if 

ar-^  ,)#2  »» 

The  general  solution  of  this  pfurtial  differential  equation  may  be  obtained  by 
noting  that  if  ve  set  0  =  —  ^  , 
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The  left  lid*  of  Equation  (12), la  than  equal  to  7 ao  that  the  aoalar  vave  Eq.  12 

^  dx^ 

1  1  iff. 


becomes 


^  ^  rC' 


The  general  solution  of  vhieh  for  outgoing  dieturbaneea  is 
^  =  F(r  -  Cpt  ).  Here,  c^  =  ( 

is  the  F-eiave  velocity.  The  displacmnent  potential  is  nov  of  the  form 

0  =  i  F(r  -  c  t) 
r  p 

Strai(^t~forvard  computation  yields  the  following  quantities: 

u  =  =  7  F'  (r  -  c  t)  -  F{r  -  c  t) 

tfr  r  p  ^2  p 

f  =4^  =  -F"  -“F*  +%F 
Cr  a  r  r  2  3 


rr 


=  F"  (r  -  c  t)  -  ^  F'  (r  -  c  t)  +  ^  F{r  -  0  t) 


°00  -  ®00  r  ^  r^  r^ 
The  maximum  stress  difference  is 


r  r 

It  can  be  seen  that  knowledge  of  stress  at  the  boundary  r  =  a  gives  us  a  differen¬ 
tial  equation  of  the  form 

AF"  (f)  -  BF'  (?)  +  CF  (f )  =  f  (^) 

where  f  is  a  known  function  and  A,  B,  C  are  positive  constants.  From  the  form  of  the 
argument  ^  =  r  -  o^t,  it  follows  that  ^  becomes  smaller  as  t  increases,  r  remain¬ 
ing  constant.  Vhen  the  disturbance  starting  at  r  =  a  at  time  zero  readies  the  point 
r,  we  have  r  -  c  t  =  a.  For  values  of  ^  ^  a  the  function  F  and  all  its  derivatives 
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Yiniih.  Vt  my  hATS  a  Btrasa  dlsooatinuiiy  ^an  tha  disturbanoa  bagina,  but  no 
diaplacamant  diaoontinulty.  Suitabla  initial  oonditiona  for  tha  problam  ara 

P*(a)  =  P(a)  =  0 

The  apparent  negative  damping  indicated  by  the  negative  coefficient  of  the  P'  term 
in  the  differential  eqimtion  for  P  is  explainable  if  one  remembers  that  ^  vill  be¬ 
came  smaller  as  t  increases.  In  time,  the  solution  damps  out  for  a  limited  distur¬ 
bance,  as  it  should. 

If  ve  no  longer  restrict  the  pressures  on  the  spherical  cavity  to  be  radially  sym¬ 
metric,  the  problem  becomes  some^diat  more  complicated.  It  is  possible,  hovever, 
to  solve  this  dynamic  elasticity  problem  in  such  a  manner  that,  as  before,  only 
the  solution  of  ordinary  differential  equations  is  required  in  the  numerical  evalua¬ 
tion.  This  procedure  avoids  the  difficult  and  tedioiis  computation  of  values  of  a 
Fourier  integral,  the  unfortunate  feature  of  the  classical  approach.  Of  course, 
in  the  problem  as  given  the  stresses  on  the  "top  half"  of  the  spherical  cavity  are 
at  our  disposal.  If  we  choose  the  stresses  to  be  antisymmetric  with  respect  to  the 
surface  (the  plane  z  =  0),  this  plane  vill  be  loaded  by  only  shearing  tractions. 

If,  on  the  other  hand,  ve  choose  a  symmetric  distribution,  the  plane  is  loaded  by 
only  normal  tractions.  In  either  CMe,  these  loads  must  be  removed  by  other 
methods. 

The  derivation  vill  be  given  in  general  terms  and  then  specialized  for  tvo  cases. 

The  first  of  these  vill  be  the  radially  symmetric  case  already  treated;  the  second 
vill  be  the  lowest  mode  of  the  antisymmetric  case.  In  vfaat  follows,  Cartesian 
Tensor  notation  vill  be  \ised  throuj^out  the  computations.  A  sum  over  n  vill  be 
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meant  only  if  the  sum  sign  occurs  explicitly. 


If  is  the  displacement  vector,  the  Navier  equations  of  elasticity  read 


/  -v  V  ^  A  W  1  /)  C/  “1 


idiere  d,  the  dilatation  is  -r — \fe  first  express  the  displacement  vector  as  the 
sum  of  tvo  vectors,  the  first  corresponding  to  no  rotation,  the  second  to  no  dila¬ 


tation. 


1  dx^ 


it/, 


■Ijk 


'  X. 
3 


tdiere  6^.  ..  is  the  alternator.  Substituting  into  the  equations  of  motion,  ve  have 
X3K 

=p[i4- (44) 


■ijk  dx,  '  ^^2 


These  three  equations  are  satisfied  if  (  7l  +  =  p0  =  P 

7(  +  2u  1/2  u  1/2 

The  constants,  ( - p  and  (^)  =  c^,  are  the  propagation  velocities  of 

F  and  S  vaves,  respectively.  Ve  see  that  the  problem  has  been  reduced  to  solving  a 
scalar  and  a  vector  vave  equation. 


Consider  first  the  scalar  vave  equation.  A  convenient  vay  of  solving  this  equation 
in  a  form  suitable  for  vise  idien  spherical  boundaries  are  involved  is  to  note  that 
the  relevant  solutions  of  the  "homogeneous"  equation  \J^0  =  0  is  0  =  K^,  lAere 
is  a  solid  spherical  harmonic  of  order  n.  Three  properties  of  the  solid  spherical 
harmonics  vill  be  used.  The  third  one  vill  not  be  needed  until  almost  the  end  of 
-the  computations. 
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1 


•  “  0  “  ■  0 


2.  !■  homogCMOua  of  order  n  in  x^,  and  so  that  by  Euler's  theorem  on 
homogeneous  funotionS|X^  =  nK°« 

3.  The  are  orthogonal  over  the  unit  sphere.  Over  a  sphere  of  radius  a  ve  have 

J'  K°K™ds  =0,  m  /  n 


I  '^‘•=STT* 


471  2&  +  2 


m  s  n 


Ve  nov  use  the  idea  of  variation  of  parameters  and  seek  a  solution  of  the  equation 

c  2 
P 

in  the  form  f  (r,  t)  E^,  vhere  r^  =  x.x..  Ve  shall  find  the  restrictions  on  the 
n  11 

functions  f  by  substituting  this  solution  into  the  vave  eqmtion  for  0.  Por  this 


purpose  ve  requiri 


^ 


1  1 


^0  ^^n  *i 
dx.  or  r 


+  f 


n  Ox^ 


^  ^  AK  *i  *i  .41  .  iii  K  "  *’• 

<^*i^  a 


n  ^  ^  k“  + -  4-^ -i -i  k“ 

2rr  rdr 


dr  ^2  r 


+  2 


^r  \^7x7  ^n  axj^a  x^ 


vhere  6^^  is  the  Kroneoker  Delta,  aero  ifi0j,lifi  =  j  (6^^  =  3).  Using  the 
properties  of  the  ve  have 
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7^0  =  k“ 


AJs.  2(1  +  n) 


The  vave  eqxiation  nov  becomes 
,  /2. 
le 


r  ^  2(1  +  n) 

1  +  „  V  _ 


(Jr 


dr 


A 

■  .  2  at^  ’ 


and  is  satisfied  provided  that 

^  2(1  +  n)  1  ^ 


d\ 


dr 


c  ^t 
P 


(13) 


Eqviation  (13)  is  a  linear  hyperbolic  partial  differential  equation  in  the  tvo  variables 
r  and  t  of  a  form  vhoae  general  solution  is  known.  For  n  =  0  we  find 

p 

the  solution  of  which  is 


f  i  p(r  .  0  t)  +  i  P  (r  +  c  t). 
or  p  r  1  p 

Ve  shall  only  retain  the  first  term,  which  represents  an  outward  moving  disturbance. 

It  may  be  verified  that  if  f  is  a  solution  of  Equation  (13), then  ~ (f  )  is  a 

n  ’  r  0  r  n 

solution  of  Equation  (13)  for  n  replaced  by  n  +  1,  It  follows  that  the  part  of  the 
general  solution  of  Equation  (13)  we  are  seeking  is 

f  (r,  t)  =  (~"^)  ~  F(r  -  c  t)  (14) 

n  ’  r  or  r  p 

We  now  turn  to  an  examination  of  the  vector  wave  equation  V  ^  As 

2  Cg 

before,  we  first  consider  solutions  of  the  equation  y  ^  ~  ^  then  multiply 
by  g^(r,  t).  The  choice  of  solution  is  governed  by  the  condition  that,  for  a 
loading  of  ihe  sphere  symmetric  about  the  X-  axis,  the  displacements  lie  in  the 
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plants  of  the  mtridlans.  As  will  be  sem  later,  a  oonTtnicnt  form  of  the  solution 


is  the  one  shown  belov: 


=  g^(r,  t) 

where  is  a  solid  spherical  harmonic  of  order  n<  Ve  prove  first  that 

-2,  ^  « 

V  ( C  .  X.  )  =  0. 

''  '  13k  1  dx/ 

4r  ‘1  =  '^idk  %  *  *i 

iV 


That  is 


and 


^2 


^  ^ijk  *i  ^Xj^  ^ijk  '  i ^ 

Furthermore,  it  is  easy  to  see  that  x^  is  homogeneous  of  degree  n  in  x^, 

3 

x^  and  x^«  It  follows  immediately  that  satisfy  equations  similar  to  Equation  (13). 

2_  V  _ 


)  =C4  4I.  (~I  '  N-  +  0. 


•)  =  0. 


^  gg  2(1  +  n)  _  1  ^  ^ 

As  before,  the  solution  of  interest  is 


8n  7 


(15) 


(16) 


For  an  axially  symmetric  problem,  this  is  all  ve  need  for  the  vector  potential.  Ve 

-  ^  ^  +  E 

u;  =  jT—  +  C  .  ...  . 


may  now  compute  the  displacement  vector 


'ijk 


Computations  will 


“A  ax.  . 

be  carried  out  for  a  single  n.  For  the  final  solution,  of  course,  a  sum  on  n  will 


be  required.  The  only  part  of  the  cooptation  that  is  not  self  explanatory  is  the 
use  of  the  identity 

6  V  =  (5-  -  6.  6.  ) 

1  jk  ^  mpk  im  3p  ip  3m 

This  last  relation  together  with  the  first  two  properties  of  the  ^  given  above 
leads  to 
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-  (n  +  1)  -  r 


Ty 


Hi*  (ii»tb«matioal)atraln  oompon«itt  may  nov  b*  oaleulatad: 


1  /  .  jK 


ij 


a- 


h‘* 


X.  X. 

i _ a 


k“  +  (n 


^g;  6.  .  X.  X. 

°  ^  .  _ n\  /  iJ.  1  .1 


7r^T7 


.)  (_11.^)k“ 


^  f-  «)8L  6^1  ^  ^  fi  i^l 

[r  Sx.  ■‘■r  ax^j 


1  I  -*n  ,  ,  ,,  ^  «a 

"  ^gL''  '2 

'n  ’  (“  +  1)  «n  -  r  — 


^  k“ 


A  =  ^ — r-^K“ 

c2  ^2 


The  stress  compohents ' are  found  frbm'the  strains  by  applying  Hooke's  Lav  for  an 

isotropic  solid: 

a.  .  =  ^  A  6.  .  +  2u£.  . 
ij  la  “  ij 

We  are  especially  interested  in  the  tractions  on  the  surface  r  =  a.  The  mit 

*i 

normal  is  jxist  — — 

r  =  a 


+ 


+ 


2 (-2^'^  + 

r 

[fn  -  (n  +  1)  g^  -  rg;] 
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Ve  St*  that  the  traction  is  in  tip  form 


^tre 


>?  =  --^’4  -  ^  K  *  0  -  -  “r,  -  <>  r 


(17) 


=-^f  5  (a;  -  2<  -  <)  *  -(.*!)«;.- 1^) 

Her*  primes  denote  differentiation  vith  respect  to  r. 

The  known  traction  on  the  boundary  is  a  normal  pressure  a  (0,  t)  in  the 

r  =  a 

axially  synnetrie  ease.  It  is  convenient  to  eiqiaad  the  stress  in  i^erieal  hasrmonies. 

(In  ^^t  follows',  the  solid  harmonics  will  correspond  to  surface  haxmonios  of  sonal 

type  only  by  virtue  of  the  axial  symmetry) .  The  expansion  reads 

a  (a,  0,  t)  =X!  •  (*)K^  on  r  =  a. 
rr  n  ^ 

Multiply  both  sides  by  and  integrate  over  the  sphere  of  radius  a, 

[  c  (a,  0,  t)  k“  ds  =  Z  *  (t)  [  kV*  ds  , 

-'s  "  “  “  Ja 

a  a 

By  the  orthogonality  properties  of  the  ve  have 


*_(t)  = 


m 


4ma 

Prom  Equation  (17)  we  have 


fll  +  1  f 
2nrf2  J 


a  K-  ds 
rr 


I  ^  -IS'  =  V  "rr  -  I  *»<*>  ^  ^  <“> 

n  j  n 

That  the  vector  equation  shoiild  yield  solutions  for  and  may  be  seen  by  noting 


that  is  a  radial  vector  and  -v —  lies  in  the  meridional  plane  if  &  corresponds 
r  o  X. 

0 
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to  zonal  harmonics.  We  are  then  dealing  vith  just  tvo  components  at  each  point  so 
that  tvo  functions  should  be  determined.  We  expect  that  =  0  since  it  multiplies 
the  only  vector  In  the  equation  not  in  the  radial  direction.  Formal  computation 
gives  us  this  result  and  assures  us  that  the  modes  decouple. 


X. 


Let  us  multiply  Equation  (18)  by  ^  k”,  sum  and  integrate  over  the  sphere  of  radius 


X.  X. 


KV'ds 


^-^K“ds 

8  <1 

a 


.T®  4t[  2m  +  2  .  ,jn 

^  (2m  +  1)  *  ^  “  (2m  +  1)  ^ 


4x  2m  +  1  /.X  471  2m  +  2 

=  ’  ‘m<*)  (am)* 


x“ 


or 

We  also  multiply  Equation  (18)  throu|^  by  ,  sum  and  integrate  over  the  unit 

0 

sphere.  For  this  purpose  ve  shall  require  the  surface  integral, 


1 


TxT  Jx 

J  V 


ds  . 


We  note  first  that  this  svcrface  integral  is  simply  related  to  the  integnl  of  the 

same  product  over  the  solid  sphere  by  virtue  of  the  homogeneity  of  the  functions 

,a 


Jv  Jn  » 


“  “j  ^*j 


xn-Hn-2  /r^2  j  f 


ds 


._s_r  f 

n+m+l  j  <)  X .  <3  X . 
^sa  3  3 


However,  the  volume  integral  can  be  evaluated  by  the  divergence  theorem: 
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a  a 


Ve  see  than  that 


Using  this  rastilt,  va  find 


4ii  2iih-1  .  ^_2m  /.  v  4to  2infl 

"T"(35r)* 


(19a) 


Equations  (19)  and  (19i)  nay  ha  solvad  simultanaously  to  give 
N?  =  0,  N?  =  -a  (t) 

«  X  w 

Thasa  rasxilts^togathar  vith  Equation  (17),proyida  us  with  tvo  ooiqplod  ordinary  dif- 
farantial  aquatioiu  if  va  racall  the  forms  of  tha  solutions  for  f  and  g,  Equations 
(14)  and  (16). 


For  n  =  0,  va  hara  K  =  1  and 

'7ir 

0  ^ 

The  Nj  equation  disappears  since  =  0. 

'o  = ;  -  %*' 

P'  p 


r  =  a 


-23- 


^0  +  ip%2_ 

>2  -  2^  +  3^ 

0  r  r  r 

This  gives  us  precisely  the  ssae  ve  derived  before. 


For 


n  =  1,  the  lowest  entisynnetric  mode,  ve  have  (or  z).  It  will  be 


remembered  that  z  =  0  is  the  plane  boundary.  Ve  have  that 

it/2 

J  j  °rr^^’  ““  ^  ^  > 

*a  ti/2  ® 


‘i(^)  =  r 
0 


a  (0,t)  sin  0  cos  0  d0 
rr 


The  solutions  for  f  and  g  are 


f  = 


_iJL 


dr 


-  F(r  -  0  t) 
r  P 


1  d  1 

8  =  7T7  7 


G(r  -  Cgt) 


The  differential  equations  expressing  the  boundary  condition  read 


rr 


r  2  3  4  2  3  4 

r  r  r  r  r  r 


=  +  aj^r  at  r  =  a 

-  f"  +  ^  f'  -  F  +  ^  g'"  -  g"  +  ^  g'  -  ^  G  =  0  at  r  =  a 
r  r  t  t 

Stresses  may  be  expressed  in  terms  of  the  strains.  The  most  interesting  stress 

components  are  the  shear  on  z  =  0,  idiich  must  be  later  removed jand  the  stress 

a--  or  a  directly  beneath  the  center  of  the  sphere. 

zs 


"31 


=  a  =%p''  .^p'  +%p  _21g"'  +  ^g" 
rtnZr2  3  4r  2 

z  =  0,  y  =  0  r  r  r  r 

,  (r=.) 

r"^  r^ 
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“m  .a  .1|,.4Eo" 

X  =  Op  y  =  0  ^  T  r  r 

.m.0,  (r  =  .) 

r"*  r 

It  is  not  difficult  to  see  hov  the  idees  of  the  shore  derivstion  may  be  extended  to 

the  case  of  spherical  layering.  In  all  but  the  oti>er  layer,  ve  diould  hare  to 

allov  incoming  as  veil  as  outgoing  vares,  that  is,  solutions  of  the  form 
1  i  n 

(--r—)  —  F(r  +  ct).  If  there  are  s  layers,  ve  hare  s  -  1  interfaces  vhere  tvo 
r  o  r  r 

traction  components  and  tvo  displacement  coinponents  must  match  (except  vhen  n  =  0). 

At  the  inner  surface,  ve  hare  tvo  equations,  as  in  the  present  problem.  Ihe  total 
number  of  (scalar)  boundary  conditions  is  then  4s  -  2.  Counting  both  incoming  and 
outgoing  vares  ve  hare  four  F's  and  G's  in  each  layer  but  the  outer,  and  tvo  in 
the  outer.  The  method  vill  then  give  4s  -  2  equations  in  4s  -  2  unknovns.  It  tdiould 
be  noted  that  these  are  coupled  ordinary  differential  equations  vith  given  initial 
values,  a  type  of  problem  veil  suited  for  solution  on  a  digital  computer. 

The  solution  to  Equation  (13)  for  a  radially  syenetric  load  vas  prograomed  for  the 
IBM  7090  computer  using  numerical  integration  of  the  displacement  potential.  Radial 
and  tangential  stresses  vere  obtained  for  tvo  media.  The  first  has  "h  =  390  ksi, 
ji  =  360  ksi  and  a  density  such  that  the  p  -  vare  velocity,  c^,  vas  6,000  fps.  The 
corresponding  Poisson's  ratio  is  0.260.  The  second  problem  has  ^  =  590  ksi, 
jx  =  260  ksi  and  the  same  density^  so  that  the  p-vare  velocity  remains  equal  to 
6,000  fps  but  V  =  0.347. 

The  radial  stress  at  various  depths  proved  to  be  independent  of  Poisson's  ratio  and 
thus  the  radial  stresses  for  both  problems  are  shovn  in  Figure  104.  The  maximum 
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t«ng«ntial  stresses  are  nearly  proportional  to  A  .  These  are  shovn  for  -various 
dep-ths  in  Figure  105  and  106  for  both  problems. 

D.  SOLUTIONS  BT  DIGITAL  SIMULATION 

Two  discrete  models  of  -the  lumped  mass-spring  type  were  used  in  the  present  study. 
The  model  used  for  the  solution  of  Problem  1  is  a  spherical  model,  vdiereas  that 
used  in  sol-ving  Problems  2,  3,  and  4  is  axially  synmetric.  The  equations  of  motion 
of  the  discrete  systems  are  identically  the  same  as  the  differential  equations  of 
the  corresponding  solid  continua.  This  is  demonstrated  in  -the  following  discussion 
for  the  spherical  model.  A  similar  demonstration  can  be  made  for  the  axially  sym¬ 
metric  use. 


In  Figure  107  is  shown  a  typical  arrangement  of  the  meuss  points  (0,  1,  2,  ...)  and 
the  stress  points  (a,  b,  c,  ...)  of  -the  spherical  model  along  the  radial  lines. 

The  masses  of  the  solid  are  concentrated  at  the  mass  points,  vdiile  the  springs  pro- 
-vlde  the  a-verage  resis-tances  of  the  solid.  The  stresses  in  the  springs  are  a-verage 
stresses  defined  at  the  Stress  points. 


Applying  Newton's  second  law  of  motion  to  a  -typical  mass  point  "0"  the  equation  of 
motion  of  such  a  mass  point  in  terms  of  stresses  can  be  deri-ved  as  follows: 


>diere: 


t  2x  ,dev‘ 


a^Nrb^) 


=  average  radial  stress  at  stress  point  "a". 

Oq  =  average  tangential  stress  at  stress  point  "c". 
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=  radial  distance  of  stress  point  "a”  from  the  center. 
P  =  mass  density 

=  radial  acceleration  of  mass  point  "0” 

J\  =  grid  distance,  see  Figure  107. 

A0  =  incremental  angle,  see  Figure  107. 


Simplifying,  the  above  equation  becomes. 


a  b  ./  o  Ox 

a  -  a  2{c  -  oj 

r  r  r  o  » 

“1 —  - 7 -  =  “o 


vhich  is  identically  the  same  as  the  central  finite  difference  analog  of  the  dif¬ 


ferential  equation. 


d  r 


2(a  -  aj 


^  ^  dt2 


This  last  eqviation  can  be  recognised  to  be  the  differential  equation  of  motion  of 
a  spherically  symmetric  solid. 


Problem  1  -  Infinite  Snace 

This  problem  involves  an  infinite  solid  space  vith  a  spherical  hole  having  a  radius 
of  660  feet  and  subjected  to  a  uniform  blast  pressure  on  the  boundary  of  the  hole. 
Results  of  the  numerical  calculations  vith  the  spherical  model  are  presented 
graphically  in  Figures  108  throu^  113.  Figure  111  shows  the  pressure-time  relation¬ 
ship  of  the  radial  stresses  for  points  along  a  radius,  idiile  in  Figure  112  are  shown 
the  same  stresses  plotted  against  the  radial  distance  for  the  specific  times  in¬ 
dicated.  Similar  pressure-time  plots  for  the  tangential  stresses  and  the  particle 
accelerations,  velocities,  and  displacements  are  presented  in  Figures  108,  109,  110, 


and  113 
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Problem  2  -  Half-Snace  with  "Correction"  Loading 

In  Figure  9Sb  is  sbovn  a  meridional  section  of  a  half-space  ^ich  is  subjected  to 
a  "oorreotion"  loading  at  a  specified  time,  t^.  This  correction  loading  is  exactly 
equal  to  the  negative  of  the  tangential  stresses  (see  Figure  113).  The  hemispherical 
cavity  is  free  of  stress  in  this  ease. 

The  stresses  and  particle  motions  derived  from  this  problem  ^en  superposed  vith 
the  corresponding  stresses  and  particle  motions  of  Problem  1  represent  the  solution 
of  the  half-space  subjected  only  to  a  uniformly  distributed  pressure  applied  at 
the  boundary  of  the  semi-spherical  cavity. 

Graphical  results  representing  time  variation  of  stresses  and  particle  velocities 
are  given  in  Figures  114  through  120.  Numerical  results  for  points  along  the  sirr- 
face  and  those  along  lines  inclined  at  15  degrees  and  30  degrees  from  the  horizontal 
are  presented.  Beyond  the  30  degree  line,  the  results  are  relatively  small  and  are 
not  presented. 

The  results  from  this  problem  are  referred  to  as  "correction''  stresses,  velocities, 
or  accelerations. 

Problem  3  -  Half-Snace  vith  Ihiform  Loading 

This  problem  involves  a  half-space  with  a  semi-spherical  cavity  with  a  radius  of 
660  feet  as  shown  in  Figure  98c.  A  uniformly  distributed  pressure  having  a  pres¬ 
sure-time  relationship  as  that  given  in  Figure  102  is  applied  on  the  entire  boun¬ 
dary  siurface  of  the  semi-spherical  cavity. 
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Plots  of  ths  radial  partiols  aooslsratlons  and  vsloeitios,  and  ilit  radial  sad  tan¬ 
gential  stresaes  aa  a  function  of  time  are  given  in  Figures  121  throu^  130i  The 
resiilts  for  points  along  lines  having  different  angles  of  inclination,  0,  with  the 
horisontal  surface  are  given)  specifically,  the  results  for  lines  vith  0  =  0, 

15,  and  30  degrees  and  for  points  along  the  vertical  line  (0  =  90  degrees)  are  pre¬ 
sented.  The  results  for  the  points  along  lines  vith  other  angles  of  inclination 
(0  >  30  degrees)  vere  not  presented  since  they  are  almost  identical  vith  those  of 
the  corresponding  points  along  the  vertical  line.  The  results  along  the  vertical 
lines  are  also  identical  vith  the  results  for  the  full-space.  This,  therefore, 
means  that  the  full  horisontal  surface  does  not  significantly  affect  the  results 
for  regions  in  the  neighborhood  of  the  vertical  axis. 

The  resvilts  of  this  problem  should  be  identical  vith  the  combined  angvers  of  Pro¬ 
blems  1  and  2.  A  close  examination  of  Figures  108  through  1-30  vill  reveal  that 
this  is  indeed  the  case  vith  the  solutions  presented  here. 

Problem  4  -  Half-Snace  vith  Non-uniform  Loading 

A  meridional  section  of  the  half-space  of  this  problem  is  shovn  in  Figure  98d. 

The  problem  is  geometrically  the  same  as  in  Problem  3  except  that  the  distribution 
of  the  applied  pressure  is  uniform  only  on  the  lover  half  of  the  spherical  surface^ 
idiile  on  the  upper  half^  the  pressure  distribution  varies  as  p^  sin  0  for  0^  0  ^  45 
degrees.  Such  a  pressure  distribution  closely  resembles  the  corresponding  results 
of  Brode  and  Bjork. 
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Results  presented  ere  for  rsdial  and  tangential  stresses,  and  for  particle  displace¬ 
ments,  velocities  and  accelerations  as  shovn  in  Figures  134  through  133.  Other 
stresses  and  particle  motions  vere  computed  but  are  not  presented  here.  The  re¬ 
sults  are  plotted  as  functions  of  time  for  points  of  increasing  radial  distances 
along  lines  vith  four  different  angles  of  inclination  from  the  horizontal. 

The  results  along  the  lines  vith  =  45  degrees  and  0  =  90  degrees  are  very  close 
to  each  other.  From  this,  it  follovs  that  the  results  in  the  region  between  the 
45  degree  line  and  the  vertical  should  be  almost  identical  with  the  corresponding 
results  on  the  vertical  line.  It  should  be  pointed  out  that  the  restilts  on  the 
vertical  line  are  also  identical  vith  those  of  the  full-space,  thus  indicating 
that  the  horizontal  boundary  of  the  half-space  has  negligible  influence  on  the 
groimd  motions  and  pressures  in  the  region  close  to  the  vertical 


E.  SIGNIFICANCE  OF  RESULTS  JWD  CONCLUSIONS 
1.  Discxission  of  Results 

The  neglect  of  a  plastic  region  between  the  hydrodynamic  model  and  the  linearly 
elastic  solid  has  an  insignificant  effect  on  the  press\ire-time  relationship  assumed 
in  this  study.  However,  since  the  study  was  predicated  on  the  assumption  of  on 
elastic  solid,  and  because  of  the  extremely  hi^  stresses  experienced  by  the  solid, 
>diich  would  invariably  cause  plastic  flow  in  the  regions  under  consideration,  the 
results  can  not  be  strictly  applicable  to  materials  idiich  are  inherently  non¬ 
linear  or  ^ich  exhibit  elastic  properties  only  at  low  pressures.  The  calculated 
ground  motions  and  stresses,  therefore,  must  be  considered  in  this  lij^t  in  relation 
to  actual  earth  materials. 
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Comparing  tho  results  of  the  full*>spaoe  vith  the  corresponding  results  of  the  half- 
space  solutions  indicates  that  the  horisontal  surface  boundary  has  negligible  effect 
on  the  conditions  in  the  lower  regions  of  the  half-space.  This  is  clearly  supported 
by  the  almost  identical  results  for  points  on  a  diagonal  line  with  those  on  the 
vertical  linOf  and  also  by  the  closeness  of  all  of  these  relationships  vith  those 
of  the  full  space.  This  fact  is  also  verified  by  the  results  of  the  "correction'' 
solutions  ^diich  show  that  the  "correction"  stresses  and  velocities  decay  rapidly 
as  a  function  of  depth.  In  the  regions  close  to  the  horisontal  boundazy,  the  tan¬ 
gential  stresses  are  significantly  affected  by  the  boundary;  the  peak  tangential 
stresses  is  increased  as  can  be  seen  fiom  Figure  127  throuf^  130,  with  its  maxi¬ 
mum  value  aome>diere  along  the  line  vith  a  15  degree  angle  of  inclination  from  the 
surface.  A  surface  effect  gave  rise  to  the  tangential  stresses  in  the  vicinity  of 
the  surface,  as  shovn  in  Figure  127.  The  radial  stresses,  on  the  other  hand,  are 
affected  to  a  much  lesser  degree  by  the  surface  boundary;  the  magnitude  of  the 
radial  stresses  influenced  by  the  boundary  are  shown  in  Figure  116  and  117. 

Since  an  expanding  air  blast  over  the  surface  vill  induce  normal  pressures,  which 
are  essentially  similar  to  the  "correction"  loading  used  in  Problem  2  althou|^ 
the  air  blast  vill  have  a  different  pressure-time  relation  and  hij^er  peak  values 
than  those  of  the  correction  loading,  the  effects  of  an  air  blast  can  be  expected 
also  to  be  similar  to  the  effects  of  the  correction  loading  and  would  be  pronovmeed 
only  in  the  regions  between  the  surface  and  the  diagonal.  In  the  lower  regions 
of  a  half-space,  the  directly* transmitted  effects  vill  predominate. 
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In  all  cases,  the  radial  relooity-time  eurres  are  similar  to  the  corresponding  radid. 
stress-time  ovarves;  hoverer,  these  tvo  functions  are  not  singly  related  as  in  the 
simple  one-dimensional  ease.  Both  the  peak  ralues  as  veil  as  the  rise  times  of 
the  radial  stresses,  Teloeities,  and  accelerations  decay  vith  increasing  radial 
distance.  The  peak  'values  of  -the  -tangential  stresses  also  decreases  as  a  function 
of  the  radial  distance,  in  addition  to  its  relations  vith  the  angle  of  inclination 
vith  the  horizontal  surface. 

2.  Decay  of  Peak  Radial  Pressures.  Velocities,  and  Accelerations 
The  rate  of  decay  of  the  maximum  radial  pressures  at  increasing  radial  distances 
are  summarized  in  Figure  133  and  156.  These  shov  that  along  -the  vertical  axis  of 
the  half-spaces,  the  radial  pressures  decay  as  the  inverse  pover  of  about  1.15  of  tie 
radial  distance,  vhich  is  precisely  at  the  same  rate  of  decay  as  the  radial  stresses 
in  the  full-space.  At  the  surface,  the  radial  pressures  decay  more  rapidly  and  is 
approximately  -vith  an  inverse  pover  of  1.40  of  -the  radial  distance. 

The  decay  of  the  peak  radial  velocities  as  a  function  of  the  radial  distance  fol- 
lovs  the  same  trend  as  the  radial  stresses  along  lines  of  the  same  angles  of  in¬ 
clination  vith  the  horizontal  siu-face.  These  are  demonstrated  in  Figures  132  and 
155.  Significantly,  these  also  shov  that  along  -the  vertical  axis  of  the  half¬ 
space,  the  peak  velocities  decay  vith  -the  same  inverse  pover  of  1.15  as  the  de¬ 
cay  of  the  peak  velocities  in  the  full-space  solid. 

The  maximum  accelerations  decrease  'vith  varying  in-verse  povers  of  the  dis-tance, 
approaching  inverse  povers  of  1.80  at  the  surface  and  1.50  along  the  -vertical  axis 
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at  graatar  radial  diatanoaa.  At  olosar  diatanoaa,  tha  daoay  taada  to  ba  lata 
rapid,  as  shovn  in  Figures  131  and  154« 

It  must  be  pointed  out  that  the  rate  of  decay  of  the  peak  values  of  the  stresses, 
velocities,  and  accelerations  for  regions  between  the  vertical  axis  and  -the  dia¬ 
gonal  is  about  the  same  as  the  rate  of  decay  along  the  vertical  axis. 

These  decay  rates  are  for  the  elastic  condition  assumed  in  the  analysis,  and  may 
be  slower  than  actual  decay  rates  in  real  earth  materials. 

3.  Specific  Conclusions 

The  following  specific  conclusions  may  be  derived  from  the  numerical  results  pre¬ 
sented  here: 

(1) .  The  region  below  the  diagonal  line  for  the  half-space  remains  essentially 

the  same  as  a  full-space;  that  is,  the  horiscntal  sxurface  boundary  has 
practically  no  effect  on  the  conditions  within  this  region. 

(2) ,  Primary  effects  of  the  horisontal  boundary  are  limited  to  the  region 

above  the  diagonal  line.  These  effects  axe  most  pronounced  on  the  sur¬ 
face  and  decrease  very  rapidly  at  greater  depths.  Tangential  stresses 
in  the  shallower  region  are  significantly  affected  by  the  boundary, 
while  the  radial  stresses  are  affected  to  a  lesser  extent. 

(3) .  The  effects  <£  normal  air  blast  pressures  can  be  expected  also  to  be 

limited  to  the  regions  above  the  diagonal. 

(4) .  The  decay  of  the  radial  presstxres  with  increasing  radial  distance  is 

less  rapid  for  an  elastic  solid  idien  compared  Irith  the  decay  in  the 
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hydrodynamie  modal  of  Brode  aad  Bjork.  Inverse  powers  of  1.40  to  1.15 
of  the  radial  distance  (depending  on  the  depth)  were  determined  in  the 
present  study  os  con  pared  with  inverse  powers  of  3  to  1.50  for  the  hy¬ 
drodynamic  case. 

(5) .  Maximum  tangential  stresses  decrease  faster  with  depth  from  the  horizontal 

surface  than  with  the  radial  distance. 

(6)  .  The  rates  of  decay  of  the  radial  velocities  are  the  same  as  the  decay 

of  the  radial  pressures.  Radial  velocity-time  curves  are  similar  to  the 
pressure-time  cvirves;  however,  there  is  no  sinple  relationship  between 
the  velocities  and  pressures. 

(7) .  The  peak  values  of  the  radial  accelerations  tend  to  decay  more  rapidly 

at  greater  distances,  with  inverse  powers  approaching  1.80  of  the  radial 
distance  at  the  surface  and  1.50  along  the  vertical  axis. 

(8) .  Rise  times  of  the  radial  pressure  pulse  decays  with  increasing  radial 

distance)  however,  no  evaluation  of  -Qiis  decay  has  yet  been  made  in  this 
study. 

F.  SraiCIAL  CASES 

1.  Non-Homogeneitv  and  Stratifications 

Certain  problems  involving  solid  spaces  with  a  number  of  layers,  each  possessing 
different  elastic  constants  can  be  treated  with  the  technique  of  digital  simulation 
used  in  the  present  study  with  sligdit  modifications  of  the  present  computer  pro¬ 
gram.  The  layeringyhowever,  must  be  restricted  by  the  following  geometric  consi¬ 
derations  : 
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(1) «  For  ipborioally  oynmoirio  probloma  (in-volving  infinito  ipaoM^for  oxu^lo) 

tho  layoring  muat  naintain  tho  condition  of  ■phorieal  oynmtryi  thia 
moana  that  the  different  layera  muat  be  apherioally  oonoentric. 

(2) ,  For  aaciallyHiynmetrio  problema  (involving  half-apaoea),  the  layering 

can  be  of  tvo  eonf igurationa : 

(a)  If  the  loada  are  applied  on  a  semi-apherioal  cavity,  the  layera 
muat  alao  be  concentric  aemi-apherea  of  given  radial  thiokneaaea. 

(b)  If  the  loada  are  applied  on  the  horiaontal  aurface,  the  layera  can 
be  horiaontal  or  apherioal. 

In  addition,  the  analyaia  muat  be  baaed  on  -the  aaaumption  that  there  are  not  relative 
motions  at  the  interface  of  any  tvo  adjacent  layers. 

2,  Spherical  Problems  in  Bilinear  Solids 

At  present  the  literature  contains  a  number  of  solutions  to  the'problems  of  plane 
vave  propagation  in  bilinear  media  (e.g. ,  Rahbmatulin,  K.  A.,  "On  the  Propagation 
of  Plane  Waves  in  an  Elastic  Medium  vith  a  Nonlinear  Stress-Strain  Lav",  Uchenye 
Zapiski,  ttiiversity  of  Moscov,  Vypusk  152,  1951,  pp.  47-555  Aksiak,  R, ,  and 
Veidlinger,  P, ,  "Attenuation  of  Stress  Waves  in  Bilinear  Materials",  Journal  of 
Engineering  Mechanics  Division,  ASCE  V.  87,  No.  EM3,  June  1961;  Sutcliffe,  S., 
"Strong  Shock  Formation  in  Bilinear  Media",  to  be  published  in  the  forthcoming 
issue  of  EMD,  ASCE). 

The  study  of  the  same  problems  in  the  radially  symmetric  system  is  not  nearly 
as  complete  but  can  be  carried  out  by  similar  analyses  (e.g.,  Rakhmatulin,  K.  A., 
and  Demianov,  Tu.  A.,  "Froohnost'  pri  intensivnykh  nagrvskakh",  Fismatgis,  Moscov 


-35- 


1961).  The  added  eomplexiiy  of  the  nonlinearity  of  the  differential  equations  in  the 
radially  synmetric  ease  makes  the  subject  some^t  more  difficult,  and  generally 
makes  approximate  analytical  techniques  necessary. 

The  differential  equation  of  motion  for  dilatational  vaves  in  a  radially  symmetric 
system  vith  linear  stress-strain  relation  is  given  by 


^^u  2  /  ^u  Uv  1  - 

(20) 

^ere  u  is  the  displacement  in  the  r  -  direction,  r  and  t  are 

space  and  time  coordin- 

j - - 

ates,  and  c  =~y  — p  ^  is  the  acoustic  velocity. 

The  characteristics  of  Equation  (20)  are 

r  ±  ct  =  R 

(21a) 

dp  i  ^  +  —  (p  —  — )  dr  =  0 
^  c  r  r 

(21b) 

The  principal  difficulty  encountered  in  the  analysis  by  characteristics  lies  in 
the  fact  that  Equation  (21b)  is  not  integrable  in  a  simple  form  unless  the  potential 
0,  vith  u  =  Is  introduced.  Thus,  as  vill  be  done  in  this  treatment,  it  is 
sometimes  more  convenient  to  use  the  general  solution  to  Equation  (20)  and  a  poly¬ 
nomial  representation. 

Let  the  stress-strain  relation  for  a  radially  symmetric  system  be  given  as  in  Figure 
157a  udiere  loading  takes  place  along  the  line  vith  slope  c'  and  unloading  along 
the  line  vith  slope  c.  Let  the  radial  stress  at  r  =  r^  be  given  by  a  function  of  tla 
type  shewn  in  Figure  157b.  In  the  (r,t)  plane,  the  solution  can  be  represented  as 
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i 

I  : 

shovn  in  Figure  157o.  Up  to  r  =  r*,  »  strong  shook  is  propegnted  along  the  line 
r  -  e't  =  0,  and  unloading  takes  place  in  region  I.  At  r  =  r*,  the  discontinuity 
.  vanishes,  and  the  elastic  loading  region.  II  forms,  vith  unloading  still  in  region  I, 


( 


1 . 


I 


Let  us  consider  only  that  portion  of  the  (r,t)  plane  for  r  <  r^,  and  make  a  series 
representation  of  the  solution  in  that  domain. 

Introducing  the  potential  0,  vith  u  =  into  Equation  (20)  yields  the  general 
solution  in  the  unloading  region 

0  =  ^  f  (r  -  ct)  +  ^  g(r  +  ct) 

u  =  i  (f  +  g')  (f  +  g) 

(22) 

-if  =  i(f*  X,«)  -^(f  +8')  +  ,) 

=  -.•)=■%(»'  -8') 

r 

Along  the  shock  front  r  -  c't  =  0,  impulse-momentvnn  yields 

c '  du  +  A  u .  =  0  (23a) 

r  V 

vhere  A  u^  and  Au^  are  the  discontinuities  in  strain  and  particle  velocity.  Or^ 
since  and  u^  are  aero  belov  r  -  c’t  =  0, Equation  (23a)  may  be  vritten 

c'u^  +  u^  =  0  (23b) 


On  the  t-axis  let  =  <?  (r^,  t)  be  given  as  a  function  of  the  type  shovn  in 
Figure  157b.  Then  letting  f(r^  -  ct)  =  f^  and  g(r^  -  ct)  =  g^,  ve  have 


a(ro, 


t) .  <»;  *  8')  -  (<  *  8;)  *  (t,  8,) 


(24) 
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On  the  shock  front  r  -  o't  =  0  Eg^unilon  (23b)  yields 


c'[i  (f*  +  g")  (f'  +  g')  +3  (f  +  g)]  + 

c[7  (f"  -  g")  (f'  -  g')J  =0  (25) 

r 

and,  since  u  =  0 

7  (f'  +  g')  (f  +  g)  =  0  (26) 

r 

Equation  (25)  becomes: 


Furthermore,  at  r  ,  0,  u  =  0  and  a  is  discontinuous,  thus 
o  r 


i  f:(r^)  -^f(r  )  +^g'(r  )  -■3~g(r  )  =  0  (28) 

r  o  ^or  0  ^0 

or  or 

0 

Equations  (24),  (27),  and  (28)  are  sufficient  for  the  solution  in  region  I  for 
r  <  r*. 


Let  t),  f  and  g  be  represented  by  the  polynomials  as  follovs: 


2  3 

a(r^,  t)  =  p^  +  pjt  +  P2t  +  pjt  +  - 

p^t 

t(r^  -  rt)  =  f(f )  -  *  f  *  - 

n^ 

g(r^  +  ot)  =  g  ( <1  )  =  gg  +  gi^  +  82^  ^  - 

f  Frd  -  ^)J  -  f(a  r)  =  f^  +  fj^a  r  +  f2a^  r^  + 

c 

.  n  n 
— —  far 
n 

g  l^rd  +  ^)j  =  g(pr)  =  gg  +  g^Pr  +  g2P^r^  +  — 

-  ggp  r 
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Substitution  in  Equation  (27)  and  squating  liks  povsrs  of  r  yields: 

*1  t  ^1 


2-  p2-2 


C ' 


«n 


rc  +  1 

'  (f)“  f  n  >  3 

P  “ 


Lc.  _S^-^ 


2 

n  -  n 


Equation  (28)  yields 

-  Tf,  +  2f  r  +  3f,r  ^  +  4f^r  ^  +  Sf.r  '*] 
r^l  CO  3o  4o  5oJ 

0 

-  [f  +  f.r  +  f  r  ^  +  f,r  ^  +  f  .r  +  f.r  ^ 

^2Lo  lo  CO  3o  4o  5o 

0 

+  r  [  gl  +  +  3g3r/  +  4g^r/  +  Sgjt/J 

o 

1  r  2  3  4  5i 

-  — I  «o  +  «l^o  Vo  Vo  Vo  +  Vo  j  =  ° 

*^0 

Equation  (24)  yields 

p  =  2  +  2f  +  6f_r  +  12f  .r  ^  +  20f_r  ^  +  - 1 

or  L  2  3  0  4  0  5o  J 

o 

-r^r  [f,  +  2f.r  +  3f,r  ^  +  4f  .r  ^  +  - ] 

2Ll  2o  3o  4o  J 

r 

0 

+ ^  ■*“  ^  ■*■  — . . —  1 

3Lolo2o  3o  J 

r 

o 

+  f  2g2  +  6gjr^  +  IZg^r^^  ^  +  — J 

ft  ^ 


(30) 


(31) 


(32) 
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] 

1 

-  ^  [  8l  +  +  383X^2  ^  ^ - j 


+  ^3  g^  +  BjTq  +  82^0^  +  g3r^^  + - ] 

T 

0 


-40- 


f,  =  r  <  ^  ^  [V]  - ^  [-*3'’] 

0  0 

+  ^  ( 71  +  ^)  (aogjO^)  -  ^  [ 84c’]  +  ^  [+  83®^  ] 

0  r  r 

o  0 

Equations (30)  throui^  (32),  sAion  aolvad  si]iiultan60ualy,vill  yield  approximate  values 
of  f  and  g  for  the  first  five  terms  ysdiieh  vill  probably  suffice  for  at  least  a  rou|^ 
solution. 


The  solution  can  be  extended  beyond  r*  by  noting  that  in  II  a^, 

a  is  continuous  across  the  unloading  vave  separating  I  and  II, 
r 

unloading  vave. 


u  =  o  and  r-o't  =  o, 
and  -77—  =  0  on  tha 


It  may  be  noted  here,  that  the  problem  could  be  solved  by  characteristics,  but  since 
the  boxmdary  conditions  are  set  on  tvo  lines;  t  =  0,  r  -  c't  =  0,  it  would  be 
necessary  to  set  up  a  finite  grid  in  the  (r,  t)  plane,  ^ich  vould  necessitate  Ihe 
solution  of  an  increasingly  larger  system  of  simultaneoxis  equations  for  each  suc¬ 
cessive  time  step. 

G.  RECOMMENDATIONS  FOR  FURTHER  STUDIES 

The  present  study  serves  to  indicate  that  solutions  of  a  number  of  difficult 
problems  involving  ground  shocks  are  possible  vith  the  technique  of  digital  simu¬ 
lation  employing  the  discrete  models  vhich  vere  used  in  performing  the  numerical 
solution  of  the  pree|ent  problems.  These  models  were  developed  vith  consideration 
for  handling  problems  of  continuum  dyiuuaics  involving  plastic  yielding  and  irrever¬ 
sible  flov,  as  veil  as  viscous  effects.  Hovever,  the  computer  programs  ciirrently 
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•iTftilable  are  reetrloted  to  linearly  elaetio  aolida,  althou|^  the  inolueion  of  tIs- 
oosity  for  eolving  problema  of  Voigt  eolide  does  not  involve  major  difficulty.  The 
handling  of  nonlinear  oharaoteristics  vill  require  further  research. 

Further  studies  that  can  be  performed  with  only  minor  or  no  modifications  of  pre¬ 
sent  cooqiuter  programs  are  the  following: 

1.  More  extensive  determination  of  the  stresses  and  ground  motions,  including 
shear  stresses,  circumferential  stresses,  principal  stresses,  and  direc¬ 
tions  of  principal  planes,  and  tangential  accelerations,  velocities,  and 
displacements  for  elastic  half-space  solids. 

2.  Determination  of  the  effects  of  an  expanding  air  blast  on  the  horizontal 
surface,  or  the  determination  of  the  combined  effects  of  the  direct  in¬ 
duced  wave  and  the  air  blast  wave. 

3.  Effects  of  elastic  layering.  The  layering  must  be  restricted  to  those 
discussed  earlier. 

4.  Effects  of  linear  viscosity. 

3.  Numerical  solution  of  spherical  blast  problems  in  solids  with  similar 
behavior. 
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Fi(5.  bscl  ;  Full  Space  with  Nole 


F/a.  ssb :  Half  Space  with 

Correction  Loading  ' 


Fig,  ;  Problem  5-  Half  Space 
WITH  Uniform  Loading 


WITH  Non-  Uniform  Loading 
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CHAPTER  V 


D  I  R E  C  T L Y  - T  R A N S M  I  T  T E  D  SROUND  PHENOMENA 
E L A S  T 0 - P L A S T  I  C  ANALYSIS 

A.  INTRODUCTION 

nils  study  DM  originstsd  in  ordsr  to  dsrivs  motion  snd  stress  vsluts  for  dirsotly- 
trsnsmittsd  ground  shook  in  the  nsor-omtsr  region,  or  so-eelled  plestio  sone, 
idiere  neither  hydrodynemio  nor  elestio  theories  is  speoifieally  applioable  to 
obtain  useful  estimates  of  these  laluss.  region  of  interest  extends  about 
the  crater,  from  the  ground  surface  to  the  rertioal  axis  below.  Very  close-in, 
within  the  crater,  shearing  stresses  are  negligible  in  comparison  with  the  extremely 
hi|^  pressures.  Someidiat  farther  out,  in  the  so-called  plastic  sone,  shearixig 
stresses  become  important.  Still  farther  out,  strains  become  small  enou|^  for 
elastic  or  yisoo-elastic  theory  to  be  used. 

Of  greater  importance  than  determining  the  distribution  of  stress  and  motion  throu|^- 
out  the  plastic  sone  are  the  computations  of  stress  and  energy  attenuation  as  a  re¬ 
sult  of  waTe  propagation  throu^  the  plastic  sone.  This  information  is  required 
because  its  availability,  as  inputs  for  the  elastic  models,  inmediately  gives 
greater  realism  to  the  estimate  of  motions  and  stresses  in  the  elastic  sons. 

In  order  to  satisfy  these  requirements  for  energy  attenuation  and  motion  solutions, 
it  is,  of  course,  essential  to  formulate  the  problem  as  an  elasto-plastio,  two- 
dimensional  model.  The  problem  is  a  taxing  technical  operation  involving  the 
solution  of  non-linear,  two-dimensional,  partial  differential  equations  and  re¬ 
quiring  the  use  of  complex  and  advanced  conqutational  techniques. 


-103. 


k  firit  sttp  toward  thoao  objoetivos  ha*  b**n  mad*  in  this  atudy.  A  formulatian 
of  th*  phyaioal  modal  baa  b**n  «cpr*aa*d  in  torma  idiioh  pormit  a  traotabla 
eomputational  proeadura.  Thia  porooadtira  baa  baan  taatad  and  dabuggad  aftar  tha 
uaual  amount  of  fala*  atarta  and  try-outa.  At  thia  ataga  only  a  pilot  problam 
haa  baan  eoliplatad  albait  with  a  amallar  numbar  of  dagraaa-of-fraadom  than  the 
formulation  ia  oapabla  of*  The  raaulta  demons  orate  tha  adequacy  of  tha  oompu- 
tational  approach  and  ceding  detail*  in  handling  a  raaliatic  phyaical  aituation. 
They  alao  indicate  a  vary  rapid  attenuation  with  diatanea  of  both  energy  and  atreaa. 
Modification  of  th*  ccnputational  procedure  haa  been  atudiad  and  planned  to  the 
and  of  handling  a  larger  ntanbar  of  dagr**a>of-fra*dom  and  of  apeading  th*  routine* 
Thia  ia  needed  in  order  to  obtain  better  dataila  in  th*  atreaa  and  deformation 
pattern. 

In  addition,  th*  ain^llfled  method  uaed  in  the  pilot  problem  prevent*  the  com¬ 
putation  of  the  reaidual  atraina  ^ftiioh  determine  the  preoiae  boundary  between  the 
plaatic  and  the  elaatio  aonea.  Modifications  of  the  procedure  again  ia  necessary 
to  get  a  measure  of  these  atraina. 

The  atudy  has  been  performed  by  Sr.  J.  A.  Brooks,  idio  has  been  assisted  by  V.  L. 
Frank,  Director,  Systems  Analysis,  Informatics,  Inc.,  in  the  computational  aspects 
of  the  work. 
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B.  STAIEMBNT  OF  FRQBIJEM 


The  apeeifie  objaotiva  of  thia  atiid^  ia  to  oomputa  the  reaponae  to  direotly-trana- 
mittad  grotmd  ahook  of  an  alaato-plaatlo,  aani-infinita  half-apaoa  aubjaotad  to  a 
aurfaee  a3q>loaion  of  a  nultl-magaton  vaapon.  Cloaa  to  the  burat,  praaauraa  in  the 
medium  are  axtramely  hi|^  and  the  ahaaring  atraaaaa  ihieh  the  medium  can  aiiatain 
are  inaignifioant  oomparad  vith  the  normal  atraaaaa  and  can  be  naglaetad.  Somaidiat 
fajrbhar  out,  ahaaring  atraaaaa  baooma  important  and  the  full  atraaa  tanaor  muat  be 
eonaidarad.  The  material  ia  atrained  vail  beyond  ita  alaatie  limit  and  alaatio  theory 
ean  be  uaad  only  aa  a  firat  approximation.  Still  farther  out,  alaatie  theory  or 
yiaoo-alaatic  theory  ia  quite  appropriate. 

The  region  of  interaat,  aa  shown  in  the  following  sketch,  is  bounded  by  an  inner 
hemispherical  surface  of  radius  r^,  the  ground  surface  0  =  V2)aDd  the  oxiter 
spherical  surface,  r^. 
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Th*  out«r  radiuit  •d'viao*s  with  tin*  at  »  T*lMltyi  r,  ■uah  that  it  ■•ptimtM 
th«  diaturbvd  tad  undliturbtd  r«gioiui* 

Solution  objoetlToa  art  to  dtfint  tht  wctont  of  tht  plattie  ton*  to  determino  th# 
attonuation  of  m»tgy  throu^  the  plaatie  lont,  and  to  dotoribo  ttrtttea,  ttraint, 
ditplaoamtntt  and  Ttlooitlet  vlthin  tba  plaatio  tent. 


C.  FORMULATION 


Tht  problwn  !■  non-linMur  ukL  tvo  dioMisional  and  iberafor*  raquiraa  nirnwrie*!  trMt- 
■Mat.  Ih«  tquatlooa  ara  giran  in  tarma  of  ganaralisad  eoordiuataa  rathar  thao  finita 
diffaranoa  approxifflatioiu  to  fiald  aquationa  in  ordar  to  avoid  apaad  limitati<m8  iia- 
poaad  by  apaca^tima  maah  ratio  raquiramanta  naoaaaary  for  oonvarganoa.  Tha  ralativa 
apaad  of  tha  two  mathoda  ia  not  known  and  maah  ratio  raquiramanta  wara  not  datarmlnad. 

It  waa  fait  that  tha  finita  diffaranoa  approaoh  had  a  hif^ar  ohanea  of  oon^lata  fail- 
\ira  by  virtua  of  apaad  and,  tharafora,  ahould  not  ba  tha  firat  approaoh  to  ba  attanptad. 
by  virtua  of  tha  ganaraliaad  ooordinata  approaoh,  ahook  diaeontinuitiaa  ara  smoothad 
ovar«  Alao,  tha  dagraa  of  approximation  to  tha  axaet  aolution  for  tha  fiald  aqua- 
tiona  dapanda  upon  tha  number  of  gtnaraliaed  eoordinataa  and  tha  ohoiea  of  baaia 
funetiona. 


Baaia  funetiona  and  aa  given  by  Equations  (H>5)  and  (H-6),*^  are  defined  over  the 
region  r  <  r  ^  r^ (t),  0 <  0  <  n/2  and  0<  0  <-271,  where  r,  0  and  0  ara  polar  ajiiarioal 
coordinates.  Tha  radius  is  r|  0  is  the  polar  angle  and  0  is  tha  meridian  angle.  The 
baaia  functions  ara  independent  of  0  and,  consequently,  the  equations  ara  restricted 
to  the  axially  symmetric  case.  The  inner  radius,  r^,  is  fixed  (i.a.,  constant 

for  any  given  problem).  The  outer  radius,  r^(t),  moves  outward  with  a  constant  velo¬ 
city,  V.  Oisplacaments  are  expressed  as  linear  combinations  of  the  basis  functions 
as  given  by 


(1) 


^Latter  H,  idiara  it  appears  in  equation  numbers,  denotes  Appendix  H. 
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(2) 


vher*  eoeffioienti,  bcoont  th«  gmcralised  ooordia*t««. 


Substitution  of  Equations  (1)  sad  (2)  into  tbs  familiar  sxprssaion  for  kinatie 
snsrqy, 


yields: 


7"”  fi^  f 


(3) 


(4) 


m  n 


idlers  A_,  B _ ,  and  F _  are  time  dependent  eoeffieients  and  P  is  the  mass  density 

on  OBI  on  \ 

uoeffioi 

/ir‘^ 


of  the  medium.  Coefficients  A  are  oi-ven  by: 

on 


(5) 


By  introducing  the  variable,  given  by: 


r= 

'I  'a 


(6) 


and  by  requiring  that  and  be  functions  of  0  and  only,  coefficients  A^  take 
on  the  form  given  by: 


o  ''o 

Tf' 

oo 


(7) 
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TIm  initgrala  in  Bquntion  (7)  nr*  ind*p*nd*nt  of  tin*  uid  oaa  b*  *T*lunt*d  ono*  for 
*Mh  problom.  Similar  formula*  for  *a&  ar*  giron  by  Equation*  (H-19)  and 
(H-20)  of  th*  AnMndix. 


DiBpl*o*m*nta  ar*  partition*d  into  elaotio  and  yi*ld*d  eonrponant*  vith  prim**  d*noting 
*la*tic  coiitpon*nta  and  doubla-prim**  danoting  yieldad  oooponanta,  Th*  total  dia- 
plaeamanty  or  unprinad  quantity,  by  definition  of  notition  aquala  th*  aum  of  the 
primed  and  double-primed  diaplaoementa.  Ih*  unprimad,  primed  and  double-primed 
notation  ia  alao  extended  to  generalized  ooordinatea  and  atraina. 


StreaaaB  are  function*  of  elaatie  atraina  given  by  th*  uaual  elaatie  relation*  be- 
teeen  atrea*  and  atrain,  Conaequently,  the  elaatie  atrain  energy  of  th*  medium  ia 
a  quadratic  function  of  the  elaatie  componenta  of  the  generalized  ooordinatea  aa 
given  by:  W"  «  /  / 

m  h 

Coeffieienta  D  ar*  funetiona  of  time  given  by  Equation*  (H-27)  throiuh  (H-30)  of 

SKI 

the  Appendix. 


The  potential  of  the  external  force*  ia  given  by; 

% 

U/-  a^(o,d)  J(t> 

O  0 

o  ^ 


(9) 
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idler*  F(t)  Uid  p(t)  KT*  prtitvr**  on  th*  •phorloal  surfae*,  r  »  r^,  aad  on  ih*  ground 
■urf»o*  r*ap*OtlT*ly«  Th*  pr***ur*  funotioui  F(i)  and  p(t)  ar*  ainq^lifiad  to 
funotiona  of  t  alon*  by  rlrtu*  of  th*  apeoifio  problema  to  b*  vorksd.  j^ydrodyxuunio 
emqnitationa*  indioat*  that  at  eartain  auitabl*  value*  of  initial  time,  t^,  moat  of 
th*  energy  tranamiaaion  aeroaa  a  suitably  ohoaen  inner  radiua,  vill  have  already 
occurred  and  practically  no  energy  tranamiaaion  vill  have  occurred  across  a  spherical 
surface  of  radius  tvice  that  of  r^.  At  the  same  time,  pressures  p(r  ,  t)  and 
F(0,  t)  vill  have  decreased  to  more  than  an  order  of  magnitude  less  than  those  mid- 
vay  betveen  r^  and  2r^.  Ih.*  air  vave  front  vill  have  traveled  ahead  of  th*  radius, 
r^,  and  the  pressure  distribution  over  the  ground  surface  from  r^  to  r^  vill  be 
approximately  constant  at  any  given  time.  Ihe  pressure  over  the  inner  surface  vill 
have  dropped  to  the  same  order  of  magnitude  as  the  ground  surface  pressure  and,  by 
virtue  of  its  insignificance,  can  be  assuMd  to  be  uniformly  distributed. 


By  making  suitable  substitutions  of  expressions  already  developed,  Equation  (9) 
becanes ;  _ _ _ 


(10) 


M 


where 


-4^ 


-  y„  fr, ?) y  Js  <“) 


*Brode,  H.  L.  and  Bjork,  R.  L. ,  "Cratering  from  a  Megaton  Surface  Burst",  RM-2600, 
1960,  The  RAND  Corporation. 
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From  tho  famlllmr  form  of  tho  Lagroagion  oquationo  of  motloa,  or  from  tho  oquiTalaat 
Hamlltoa'a  prinoiplo,  tho  oquatlono  of  motion  art  foimd  to  bat 


t[  n  n 

irtioro  oooffioionta  A  ,  B  ,  C 
mu  mi  mi 

olaotio  oaaOf  in  lAiiob 
for  dotormination  of  tho  q,"  , 


Wis  IjZ. . . . 

ai]d  aro  givon  in  Appondix  H.  Exeopt 
an  auxiliary  oot  of  yiold  oquationo  must 


(12) 
in  tho 

bo  doYolopod 


Suppose  wo  apply  Coulaod>'8  equation,  which  is  given  by  Equation  (13),  and  dotormino 

^// 

values  of  yield  strains*,  “d  c^,  such  that,  for  given  relations, 

the  following  expression  is  satisfied: 

/2'/^  (13) 

idiero  t  is  the  shearing  stress,  e  is  odhesion,  a  is  the  normal  stress  and  0  is  the 

angle  of  internal  friction.  Lot  us  denote  the  values  of  and 

rr  00’  00  r0 

idiich  have  been  determined  from  application  of  Equation  (13)  by  f^,  f^,  f^  and  f^ 
respectively.  Then  let  us  seek  values  of  q^^^,  n  =  1,  2  ....  such  that  yield  strains 
computed  from  the  q^'^  approximate  the  values  given  by  f^,  f^,  f^  and  f^.  Fy  substi¬ 
tution  of  yield  displacements  determined  from  formulas  such  as  Equations  (l)  and  (2) 
into  formulas  for  strain  we  obtain: 


(14) 


*The  strain  notation  is  fairly  obvious.  Only  the  expression  for  f  .  mi^t  cause 
some  confusion I  it  is 
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(15) 


n 

(16) 

'i  y- +^]ir|  (17) 


Then,  if  ve  require  the  following  equation  to  be  eatiefied  for  ar>itrary  yariations 
of  streas,  6o^,  6a^,  So^  end 


+  2(6/^  -i^)£ 9  Jlyel9  =  O  ^ 


(18) 


we  obtain  the  following  set  of  aimultaneoua  eqviationa  for  determination  of  the  q ' 


*"  ' 


(19) 


idiere  G  is  found  from: 
m 

JVi 

'0 


ft  )C -ff)  +  2^ 


(20) 
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(20  oon't.) 


The  ootffioicnia  in  Equation  (19)  ar*  the  aaiM  aa  tha  D_  eoaffieianta  in  tha 

ami  m 

aquaiiona  of  motion.  Hanea  iha  aquationa  of  motion  baooma: 


In  tha  "aimplifiad  thaory",  for  vhich  Equationa  (21)  ara  the  aquationa  of  motion^ 
atraaaaa  dapand  only  on  tha  inatantanaoua  total  atraina  and  not  on  tha  paat  atrain 
hiatory.  The  more  conq>lata  theory  ia  contained  in  Appendix  H.  In  place  of  Equa¬ 
tion  (19)  I  tha  more  conqtlata  aquation,  copied  from  Equation  (H-97),  ia  aa  followa: 


Equation  (22)  haa  bean  coded  into  tha  exiating  machine  program,  but  temporary 
inaocuraciaa  in  tha  E^  matrix  render  tha  more  conplata  coda  inoparatiya. 
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It  rwuftins  to  prorido  fonmilM  for  oainput*tion  of  f^,  f2»  f^  oad  f^.  In  doing  so 
VO  prootod  to  ikoteh  tko  dotail  eontninod  in  tho  Appondix  quito  sporooly.  Coinbining 
tho  equation  of  Mohr'a  eirelo  vith  Equation  (13),  vo  obtain: 


(23) 

(24) 

(25) 

idioro  and  are  principal  streasoB  ordered  auoh  tbat  ^  4!  c  cot  0, 

Equations  (23)  throu^^  (25)  are,  in  reality,  only  more  explicit  farms  of  Coulomb's 
equation. 

Ve  nov  compute  the  principal  strains  vhicb  for  the  sinqplified  theory  are  giyon  by: 


(26) 

(27) 

^3- 

(28) 

For  the  more  oonqplete  theory,  principal  strains  are  giren  by  Equations  (H-77) 
throu^  (H-79). 


-  114  - 


Next,  l«t 


f/=  5  ^  g  (-^J  -  §(^^) 

(29) 

(30) 

e:- 

(31) 

^sr#  eooffioients  § q,  “id  are  to  bo  dotonninod  auoh  that  Equations 

(26)  throui^  (28)  are  satisfied.  It  is  foimd  that  if 

C  Ca^  ^  Q(  4.  ^(£,  ^  ^ 

(32) 

then 

(33) 

If  Equation  (32)  is  not  trut,  then 


?.=  I- 


C 


Haring  ^q,  if  ^  6^^  and  if 


(34) 


+  Zyx( 


(35) 


than 


f,  ^  Sz  ^ 


(36) 
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Othervite,  if 


Ci  -  (^) 

and  if 

than 

_  /  _  C<^4> 


•od  ^2  =  ^3  =  '^*  conqplata  sat  of  poaaibilitiaa  includas  8  othar  casaa  vhich 

^'1  ^'1 

ara  given  in  tha  Appandix.  Having  obtainad  and  Gy  a  tranaforaaiion 

invaraa  to  Equationa  (26),  (27)  and  (28)  givaa  f^,  f^,  f^  and  f^  for  uaa  in  nunari- 

cal  avaliiation  of  G  ,  m  =  1,  2 
m’  ’ 


Tha  machine  program  vaa  originally  coded  from  vorking  reporta  idioae  contanta  are 
eaaantially  tha  aama  aa  given  in  Appendix  H  of  thia  report,  except  for  differencea  in 
notation  and  grouping  of  terma.  Coding,  formula  evolution,  and  reporting  never 
reached  exact  coincidence.  Tha  code  waa  gradually  "de-buggad"  and  modified  in 
order  to  eliminate  trouble  idiich  manifeata  itaelf  aa  a  violation  of  the  principle 
of  oonaervation  of  energy.  In  addition  to  correction  of  coding  errora,  tvo  modi- 
ficationa  had  a  prounoueed  effect  on  the  oomputationa. 


Tht  firit  modifleation  eoDO«ra«d  th*  •iflult*a*otui  int«gr»tl<m  of  tho  oquatioiui  of 
oiotlon,  girm  by  Equation  (12),  and  iht  atopwiso  adTaneanant  of  tha  yield  ooaponanta 
of  the  generaliaed  ooordinatea  by  use  of  Tquation  (22).  The  prograa  uaea  a  Runge- 
Kutta  integration  routine  idiioh  requiree  oonputation  of  derivatiyea  at  timea  t, 
t  +  h/2,  t  +  h/2  and  t  +  h  for  integration  over  the  interval  from  t  to  t  +  h. 
Aaaoeiated  vith  aueoeaaive  derivative  eooqputationa  vithin  an  integration  eyele,  tho 
values  of  At  in  Equation  (22)  were  originally  taken  to  be  0,  h/2,  0,  and  h/2, 
corresponding  to  time  intervals  from  t  to  t,  t  to  t  -t-  h/2,  t  +  h/2  to  t  h/2  and 
t  +  h/2  to  t  +  h.  The  original  sequenop  of  At  values  in  Equation  (22)  was  changed 
to  0,  h/2,  h/2  and  h,  corresponding  to  time  intervals  from  t  to  t,  t  to  t  -f  h/2, 
t  to  t  h/2  and  t  to  t  +  h.  This  ohange  produced  physically  plausible  results  for 
short  test  runs  but  energy  growth  trouble  returned  when  attenpts  were  made  to  inte¬ 
grate  over  durations  of  time  suitable  for  practical  problems. 

The  second  modification  was  made  when  it  was  learned  that  inaccuracies  in  matrix 
elements,  E^,  wore  causing  the  difficulty.  The  change  consists  of  replacing 
Equation  (22)  by  Equation  (19)  and  replacing  Equatioxu  (H-77)  throu^  (H-79)  by 
Eqxiaticns  (26)  throu^  (28).  Errors  in  the  E^  elements  are  on  the  order  of 
(At)  by  virtue  of  its  derivation.  This  leads  to  errors  in  the  time-wise  advance¬ 
ment  of  Equation  (22)  idiich  are  not  acceptable.  Ihis  inadequacy  can  be  corrected 
in  the  derivation  and  ineprporated  into  the  code,  but  lack  of  time  prevmted  doing 
so  in  this  study. 
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D.  PILOT  PROBLEM 


A  pilot  problom  h*i  boon  run  ^leh  dwoonotrotoo  tho  odoquaoy  of  tho  mothod  ond 
oorrootnooa  of  tho  eodod  proooduro.  Tho  problem  porojnotoro  are  m  follows: 


Nvmiber  of  degrees  of  freedom 

N  =  4 

Lamp's  constants 

A  =  301  ksi 

p  =  431  ksi 

Density 

^  =  3.3  slugs/ft' 

Cohesion 

c  =  100  psi 

Dilatational  wave  velocity 

7250  ft/sec 

Air  pressure  parameters  for 

2  MT  burst 

Inner  radius 

r^  =  250  ft 

Initial  outer  radius 

r, (t  )  =  300  ft 

1  0 

Beginning  of  cmoputation 

22  milliseconds 

Initial  strain  energy 

sero 

Initial  kinetic  energy 

3.07  X  10^^  ft  - : 

Initial  conditions  wore  ehoson  to  giro  a  volocity  distribution  which,  qualitativoly, 
had  tho  appoaronco  of  tho  Tolocity  distributions  shown  in  RAM)  report  number  RM-260O 
by  Erode  and  Bjork«  Tho  magnitudes  of  tho  velocities  were  chosen  such  -Uiat  the 
initial  kinetic  energy  was  100  KT  (TNT  equivalent) .  The  initial  strain  energy  was 
taken  to  be  soro  since  code  cheoks  had  shown  a  rapid  re-apportionment  of  energy 
irrespective  of  the  initial  assunptions.  The  initial  conditions  are  as  follows: 


7^3 

■ff/Sec. 

(t.)  = 

7^3 

St/ Sec. 

fsCt.)  = 

7^3 

St/ Sec. 

Si /Sec. 

Thla  ohoie*  giTts  •  p««k  initial  horiioatal  ralooity  at  tha  ground  turfaeo  of  about 
2.6  foot  par  milliaaoond  and  a  paak  initial  vartieal  Talocity  (not  at  tha  aana  point) 
at  tha  ground  aurfaea  of  about  1.2  foot  par  milliaaoond  inward.  Diraetly  balov  tha 
burat,  tha  paak  ralooity  ia  about  3.5  foot  par  milliaaoond  at  tha  baginning  of 
oomputationa . 

Raaxilta  obtainad  from  tha  eonqputar  ara  tabulatad  at  ragularly-apaead  intarvala  in 
tima,  angla  0,  and  coordinata^j  thay  ara  than  plottad  againat  radiua  for  fixed 
Taluaa  of  tima  and  angla  0.  From  thaaa  plota,  quantitiaa  can  ba  ra^lottad.a8 
funotiona  of  tima  or  angle  0* 

Figure  158  ahowa  radial  atreaa  veraua  time  for  rarioua  radii  along  a  line  idiioh 
makea  an  angle  of  72  degreea  with  the  vertical  axia.  Figure  159  ahowa  radial 
atreaa  veraua  tima  for  pointa  along  tha  vertical  axia.  Invaraion  of  atreaa  mag- 
nitudea  from  idiat  one  would  es^ct  ia  due  to  too  large  a  value  of 
initial  oonditiona.  the  initial  velocity  diatributlon  for  the  pilot  run  aiaqply  doea 
not  match  the  intended  initial  diatributlon  and,  aa  a  reatilt,  tha  ragion  balov  tha 
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or»t*r  ia  rapidly  put  into  tanalon.  Tha  Uniting  Taluo  of  tanalon,  «hieb  la  o  oot  0, 
abova  up  aa  aaro  on  tha  aeala  of  Flguraa  158  and  159>  Of  oouraa>  mora  appropriata 
initial  oondltiona  oan  ba  uaad. 

Reaulta  ahovn  in  Figure  158  and  159  vara  ao  diaturbing  that  hand  oomputationa  vara 
performed  to  apot  ohaok  tha  machine  computation  of  yield  atraina.  An  error  vaa 
found  but  it  could  not  poaaibly  hava  influenced  tha  reaulta  appreciably.  Tha  error 
had  tha  affect  of  limiting  the  tanaila  atraaaaa  to  approximately  600  pai  inataad  of 
172  pai,  aa  had  bean  intended.  Figure  160,  idiioh  ahova  time  hiatoriaa  of  tangential 
atreaa  for  0  =  72  dagreaa,  ia  given  for  conpariaon  vith  Figure  158.  Alihoiif^  tha 
time  acalea  uaed  in  Figurea  158  and  160  are  different,  it  can  be  aeen  that  the 
general  ciurve  ahapea  for  both  radial  axid  tangential  atreaaea  are  quite  almilar. 

Aa  ahovn  by  theae  figurea,  tangential  atreaaea  apparently  begin  to  attenuate  more 
rapidly  than  radial  atreaaea  beyond  a  diatance  of  1,000  feet.  No  explanation  of 
thia  phenomenon  haa  been  found  aa  yet. 

Radial  diaplacementa  along  the  vertical  axia  (0  =  0  degreea)  are  ahovn  in  Figure  161} 
radial  diaplacementa  for  0  =  72  degreea  are  ahovn  in  Figure  162.  Theae  reaulta  are 
not  a  good  repreaentation  of  the  true  phyaical  behavior.  Thia  occurrence  in  the 
computed  reaiilta  oan  be  explained  on  the  baaia  of  a  deficiency  in  the  choice  of 
baaia  functiona.  Of  the  baaia  functiona  currently  contained  in  the  computer  pro¬ 
gram,  none  haa  a  radial  atrain  component  different  from  aero  at  r  =  r^.  Thia 
deficiency  can  be  corrected  relatively  eaaily  and  aueh  a  correction  ia  an  obvioua 
neceaai-fy  before  further  cooiputationa  are  made. 
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Figure  163  ihewi  tangential  dlaplaoeaente  along  6  s  72  degrees  for  rarloua  radii* 

The  ourres  extend  to  223  mlllleeoonda.  Congputatlone  were  performed  out  to  490 
mlllleeeonda  and  the  tangential  dlaplaeemente  peak  at  approximately  280  mlllleeoonda . 
Ihe  faot  that  tangmtlal  dlsplaoemaots  are  large  Is  just  another  result  of  the 
Initial  reloelty  pattern  lAiioh  was  chosen.  Jhiother  factor  contributing  to  the  large 
dlsplaoemants  shown  in  Figures  161,  162  and  163  is  the  inability  of  the  "sinplified 
theory"  to  leare  behind  residual  displacements . 

The  oonputer  results  are  heartening  in  that  they  indicate  the  readiness  of  a  usable 
computer  code  and  in  that  they  indicate  a  hi(^  attenuation  of  stress,  HoweTer,  any 
oonolusion  to  be  drawn  from  them  must,  of  oourse,  take  into  aeoount  the  preliminary 
nature  of  the  results.  These  results  are  intended  only  as  an  indication  of  the  pos¬ 
sibilities  of  approaohing  an  extremely  difficult  and  complex  problem  in  a  greatly 
sinplified  manner.  That  certain  physical  conditions  have  been  violated  is  not  at 
all  surprising.  Correoting  such  disorepancles  is  part  of  the  usual  effort  in 
evolving  such  teohniques. 
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E.  CQNCLUSKKS  AMD  RE(X»MENMTI(»ra 

Th«  following  eomnonts  ■unanriio  tho  Brooks  study  of  the  propagation  in  soil  of  an 
axially  syanetric  wave  generated  by  a  nuclear  detonation  at  the  ground  surface. 

The  general  ob jeotiTe  of  this  analysis  is  to  provide  motion  and  stress  values  for 
the  directly-transmitted  ground  shock  in  tbe  near  crater  region  ^ere  plastic  dis¬ 
turbances  are  likely  to  occur.  These  comnents  define  the  extent  of  the  results 
obtained  thus  far  and  outline  the  work  necessary  to  carry  further  the  already- 
developed  analytical  techniques  to  the  level  and  range  of  conditions  required  for 
confident  engineering  design. 

The  model  is  that  of  an  elaato-plastie  medium  bounded  by  the  ground  surface  plane 
and  two  concentric  spherical  surfaces  of  radii,  r^  and  r^.  The  inner  surface  des¬ 
cribes  approximately  the  boundary  between  tbe  hydrodynamic  and  the  plastic  regimes. 

The  outer  surface  describes  the  bounda3ry  beyond  which  tbe  medium  is  undisturbed  at 
any  given  time.  Tbe  model  is  specifically  developed  to  treat  the  plastic  two- 
dimensional  behavior  of  the  ground  at  the  crater  periphery.  The  hydrodynamic 
inputs  are  derived  from  tbe  Brode-Bjork  data  ^icb  were  graphically  reduced  to  a 
more  suitable  and  simple  form.  For  the  conditions  computed  by  Brode  and  Bjork,  the 
radius  of  the  interior  sphere  is  approximately  70  meters  or  less. 

The  model  is  designed  to  provide  displacements,  stresses,  strains,  and,  in  particular, 
partitioning  of  strains  into  elastic  and  plastic  conqKments .  These  values  are  ob¬ 
tainable  for  any  selected  net  point  in  time  and  space  for  any  ground  medium  des¬ 
cribed  parametrically  by  the  angle  of  friction  0,  cohesion  c,  and  Lame  constants 
A  and  p..  The  net  points  are  defined  at  regular  time  intervals  At  at  regular 
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intervals  in  sphsrioal  eoordimtss  0  and  r,  shsrs  0  is  the  sngls  bstwssn 
the  vertioal  axis  end  radius  r  going  throu|^  the  burst  point. 

A  oompMting  code  is  available.  The  basic  fomulaticn  uses  a  generalised  coordinate 
approaoh  rather  than  a  finite  difference  approach.  It  incorporates  certain  sinpli- 
flcaticos.  For  instance,  shock  discontinuities  in  the  region  between  the  hemis¬ 
pherical  boimdarles  have  been  ignored  and  smoothed  over  by  means  of  a  series  ex¬ 
pansion  ^ioh  is  the  seme  over  the  i^le  region.  Ihis  siaiplification  does  not  in¬ 
validate  the  two^imensional  aspect  of  the  model. 

In  a  similar  fashion,  the  air  blast  pressure  over  the  ground  surface  at  any  given 
tine  is  assumed  not  to  extend  beyond  the  radius  of  the  outer  boundary.  This  sis^li- 
fied  assumption  does  not  allow  the  simultaneous  treatment  of  air-induced  and 
directly-transmitted  ground  shock  but  the  soil  confinement  from  the  air  load  is 
properly  simulated.  The  code  was  developed  for  the  IBM  1604  computer.  It  can 
compute  the  components  of  stress,  elastic  and  plastic  strain,  and  displacements 
at  net  points  within  the  region.  At  present,  the  code  also  incorporates  a  simplified 
yield  condition  in  idiich  the  stresses  depend  solely  on  total  strains  and  Coulomh's 
equation.  Only  one  pilot  problem  has  bear  solved  using  relatively  simple  soil 
parameters  and  weapon  inputs.  Resulting  soil  stresses  and  motions  have  been  com-' 
puted  for  only  a  few  space  points  and  time  intervals.  Because  the  small  nunher  of 
generalised  coordinates  used  in  this  pilot  problem  restricts  the  detail  of  the 
solution,  the  results  do  not  provide  final  design  data.  They  do,  however,  demon¬ 
strate  the  adequacy  of  the  computational  approaoh. 
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Th«  aaalyalB  and  eonputar  program  doTolopad  in  thia  study  art  nov  rtady  for  solu^ 
tion  of  physically  significant  surfaot  burst  probltos.  At  present,  the  numerical 
proeedure  employs  stress  ooiiq;«tatlon  routines  idiioh  do  not  depend  on  the  past 
strain  history  but  only  <m  ourrent  strains.  Thus  the  present  code  does  not  in¬ 
clude  the  ability  to  handle  residual  strains  and  displacements,  but  ihere  is  no- 
indication  -that  this  feature  cannot  be  added. 

Preliminary  results  fiam  -the  pilot  problem  used  to  check  the  programming  indicates 
that  directly  transmitted  ground  shock  from  a  2  MT  burst  vill  be  insignificant  in 
a  atructvcral  design  sense  beyond  a  radius  of  about  2,000  feet. 

Further  effort  toward  making  the  Brooks  solution  more  useful  in  the  definition  of 
imder ground  effects  logically  divides  itself  into  two  sequential  phases,  -the  first 
being  concerned  with  improving  -the  code  and  -the  second  consisting  of  parametric 
solutions . 

Improvement  of  the  code  has  tvo  concurrent  objecti-ves.  The  first  is  to  improve  -the 
available  code  by  incorporating  a  more  sophisticated  yield  condition  in  which 
stresses  depend  upon  the  past  strain  history.  The  second  is  to  enhance  the  speed 
of  computing  operations  and  thiu  substantially  diminish  the  burden  and  coat  of  the 
parametric  studies  of  soil  response  in  the  plastic  region. 

The  present  difference  eqxiation  form  of  the  yield  equations  has  been  found  less 
accurate  than  the  companion  differential  equation  form  of  the  equations  of  motion. 
This  difference  equation  can  be  replaced  by  a  hybrid  differential  difference  for¬ 
mulation  idiioh  -will  correct  -the  difficulties  inherent  in  a  generalised  coordinate 
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approMh.  Tha  ganaralisad  eoordiaata  ayatam  «aa  adoptad  aarly  in  tha  analyaia  of 
tha  modal  in  ordar  to  prorida  tha  aiaqplifioation  naeasaary  for  aolution. 

Tha  laek  of  apaad  in  tha  ooe^tational  prooadura  haa  haan  found  to  ba  dua  primarily 
to  raeonputation  of  quantitiaa  uhioh  could  bo  atorad  in  tha  cooiputar  mamory.  Thare- 
fora«  tha  ooda  ahould  ba  modifiad  in  order  to  optimiaa  data  tranafar  and  uae  of 
tapa  unita. 

Faramotrio  atudiaa  ara  naadad  both  to  indicate  -the  aanaitlTity  of  computed  reavilta 
to  obangea  in  aoil  parametara  and  to  provide  free-field  deaign  inputa  for  the 
analyaia  of  atruoturea  in  the  plaatic  aone.  Aa  a  baaic  program  to  ba  followed  it 
ia  auggaated  that  throe  typical  earth  materiala  be  ehoaen  for  atiidy:  a  hard  rock, 
a  aoft  clay  or  other  plaatic  aoil,  and  a  aoil  of  character iatioa  between  thoae  of 
rook  and  aoft  clay. 

Each  of  theaa  media  can  be  deaoribed  by  a  aet  of  parametric  valuea,  (i.e.,  ^ 

0,  end  o).  Computed  reaulta  of  atreaaea,  atraina,  diaplaoementa  and  velocitiea 
ahould  be  obtained  for  each  aet  of  parametric  valuoa  rather  than  for  independent 
mriation  of  aoil  parametara.  The  baab  inputa  to  auoh  enalyaea  are  thoae  found 
in  the  Brode-Bjork  eomputationa  of  BM-2600  (and  aubaequent  amended  and  unpubliahed 
data  uaing  a  finer  meah)aa  well  aa  any  later  oooputationa  by  AFSVC.  Theae  data 
provide  both  initial  energy  end  preaaure-veraua-time  functiona,  each  of  idiioh  ia 
atifficient  for  total  excitation  of  the  plaatic  regiona.  The  reaponae  of  one  of 
the  three  media  to  each  of  theae  functiona  ahould  be  confuted  and  compared  in 
order  to  aacertain  their  compatibility  with  the  general  mode?t.  aa  wall  aa  with  tha 
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■•t  of  parMwtrio  T»luts  dosoribing  tho  modi*. 

Thui  »  first  txtsnslon  of  tbs  Brooks  study  involTss  s  minismiB  of  four  sots  of  oooh 
putations.  It  must  be  understood  that  -tiie  use  in  the  eoinputations  of  the  RAND  data 
automatically  involves  certain  uncertainties  and  limitations.  The  influence  of  the 
plastic  sons  behavior  on  the  pressures  in  the  hydrodynamic  none  is  ignored.  This 
influence,  however,  has  been  shown  by  Newmark  to  be  small  end  the  results  of  a 
parametric  study  idll  enhance  considerably  our  quantftatlve  understanding  of  the 
physical  phencmena.  Such  a  study  should  be  primarily  aimed  toward  distinguishing 
between  the  environment  created  by  the  wave  propagation  at  the  surface  on  one  hand 
axid  at  depths  below  the  crater  on  the  other.  The  Brooks  model,  in  fact,  is  not 
desiffied  to  measure  small  looalised  variations  at  minute  depth  dlfferenoesi 
consequently,  the  rays,  (along  which  information  is  reported  as  a  fxmction  of  dis¬ 
tance  and  time)  should  be  limited  to  a  small  nunber.  The  number  of  rays  and  points 
on  each  ray  will  be  determined  by  the  reqxiirement  for  graphical  presentation  of 


the  results 
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CHAPTER  V  I 


COMBIKED  EFFECTS 


By  coiiblned  effects  Is  neant  primarily  the  combination  of  air>lnduced  and  directly- 
transmitted  ground  effects.  Another  conbined  effect,  that  of  the  Rayleigh  vave, 
is  a  result  of  coiqplex  Interactions  of  the  dllatatlonal  and  shear  waves  and  is 
all  alr-ixiduced.  It  has  been  discussed  in  Chapter  II  (Vol.  l).  Before  discussing 
means  of  confeining  air- induced  and  directly-transmitted  effects,  it  is  helpful  to 
summarize  briefly  the  state  of  the  art  of  determining  directly-transmitted  effects- 

A.  STATE  OF  THE  ART  —  DIRECTLI-TRANaaTTED  EFFECTS 
1-  Hydrodynamic  Analysis 

The  hydrodynamic  cratering  analysis  of  Brode  and  BJork  has  been  used  as  input  to 
both  of  the  analyses  of  directly- transmitted  ground  effects;  therefore,  any  un¬ 
certainties  or  approximations  in  the  hydrodynamic  analysis  will  affect  the  results 
farther  out.  Several  approximations  were  made  as  well  as  certain  assxaq^tions. 

One  of  the  approximations  is  the  equation  of  state  of  the  only  material  considered, 
Nevada  tuff;  the  approximation  is  satisfactory  for  this  material  but  not  neces¬ 
sarily  for  other  ground  materials.  However,  one  of  the  results  of  Newmailc's 
analysis  is  an  indication  that  the  computations  are  not  sensitive  to  variations 
in  the  equation  of  state  l‘'ui’ther  effort  ia  required  to  establish  this  important 
conclusion. 

Details  of  the  weapon  decoiq>osition  and  other  early-time  phenomena  have  been  based 
on  assumptions  which  are  based  on  highly-classified  weapon  detail  information; 
hence,  calling  the  detonation  a  2  NT  burst  is  open  to  question.  Knowledge  of  the 
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basla  of  these  «M\nvtlone  is  of  llttlo  v«lue  since  the  attacker's  weapons  will 
remln  an  unknown  quantity.  Bius,  a  oajor  Input  will  always  he  uncertain,  a  fact 
^ch  slsply  iBust  he  accepted.  Location  of  the  explosion  is  also  Halted  at  pres¬ 
ent  to  the  ground  sux'fhee  since  the  BASD  conputations  were  aade  only  for  a  surface 
hurst.  Depth  effects  should  he  considered  in  any  further  ooiputatlons  undertaken. 

The  Eulerlan  approach  to  the  hydrodynaalc  prohlea  carries  with  It  certain  diffi¬ 
culties  and  has  resulted  In  soas  sasarlng  of  the  shock  front.  This  Is  not  a 
critical  difficulty  xior  does  It  hrlng  ahout  posslhle  large  errors. 

In  suaaary  It  can  he  stated  that  the  hydrodynaalc  analysis  of  Brode  and  BJoxlc  Is 
highly  satisfactory  as  a  aeans  of  deterMnlng  energy  partitioning  between  air  and 
ground- shock.  Conputations  based  on  this  analysis  have  proven  to  he  adequate 
as  Inputs  tc  analyses  of  the  plastic  luid  elastic  regions  beyond  the  crater.  The 
chief  shortcoalng  of  the  commtations  lies  In  the  Halted  extent  of  the  results. 

More  aaterlals  need  to  he  studied  as  well  as  detonations  at  various  depths.  Until 
such  results  are  ohtsdned  It  la  necessary  to  Halt  the  conclusions  which  can  he 
drawn  from  the  plastic  and  elastic  analyses  of  Bdrooks  and  Mewanzic. 

2.  Elastic  Analysis  (Bewaark)* 

Nevaark's  analysis  takes  as  inputs  the  hydrodynaalc  results  at  a  radius  of  ahout  6^0 
feet.  These  are  applied  directly  to  an  elastic  half-space  without  any  consideration 
of  plastic-zone  effects,  Including  attenuation  of  peak  stress  and  change  of  pulse 
shape.  Neglect  of  the  plastic  zone  was  dictated  by  available  aeans  of  obtaining  approx- 
laatlons  of  response.  The  basic  aodel  used  by  Newaark  Is  a  aass-sprung  approxlaatlon 
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incorporating  radial  and  tangential  degrees  of  freedom.  To  this  model  were 
applied  pressure- time  inputs  based  on  the  hydrodynamic  results  An  ioportant 
result  of  one  part  of  the  Mevnax^  study  is  that  the  variation  of  material  prop¬ 
erties  from  one  side  of  the  hydrodynamic  boundary  to  the  other  (i-e.,  from 
hydrodynamic  to  elastic  in  Hevmrk's  analysis)  does  not  cause  serious  reflec¬ 
tions  or  other  difficulties  in  matching  hydrodynamic  outputs  to  elastic  inputs, 
■ecause  Irode's  results  %Mre  for  tuff,  Hevnark's  results  were  limited  to  a 
similar  material. 

Results  of  these  computations  are  extremely  Interesting  but,  as  is  to  be  ex¬ 
pected  vrlth  first  results  of  coapllcated  analyses,  raise  several  questions  that 
must  be  answered  before  practical  use  cem  be  made  of  the  results.  First,  the 
results  show  lar^e  tensile  tangential  stresses  far  beyond  the  ability  of  ground 
materials  to  i-esist.  Furthermore,  results  are  obtained  only  for  fairly  large 
depths  at  moderate  distances  from  the  burst  because  the  ray  nearest  the  surface 
emanates  from  the  burst  at  a  15  degree  angle  from  the  horizontal.  Also  signifi¬ 
cant  for  shallow-depth  effects  is  the  neglecting  of  air  blast  loading  on  the 
surface . 

In  addition  to  closing  the  gaps  in  the  results  as  listed  above,  results  from  the 
Newmark  analysis  should  be  extended  to  layered  media  aixi  viscoelastic  media. 

Ihe  analysis  is  capable  of  handling  these  more  complex  problems  with  only  mod¬ 
erate  changes  and  additions. 
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3*  Hasto-Plastlc  Analyslg  (Brooks)* 

Die  Brooks  plastic 'region  analysis  vas  undertaken  precisely  to  close  the 
betveen  the  hydrodynanle  and  elastic  regions.  Newaark's  results  extend  the 
elastic  region  to  over  100,000  pel,  ^Ich  is  far  above  the  elastic  stress¬ 
carrying  ability  of  earth  oaterlala.  The  aaln  results  ^Icb  vere  sought  from 
the  plastic  analysis  are  the  llnlta  of  the  plastic  tone  and  the  amount  of  at¬ 
tenuation  of  energy  and  stress  within  this  zone. 

This  analysis  employs  generalised  coordinates  and  numerical  approximations  and 
Includes  a  means  of  detecting  the  cessation  of  Inelastic  behavior.  Further,  It 
ellminateB  one  of  the  problems  encountered  In  Mevmark's  results,  that  of  high 
tensile  stresses,  by  limiting  tensile  stresses  to  a  level  which  a  ground  material 
can  withstand;  Couloafe's  equation  la  the  basis  of  the  assumed  behavior  of  the 
ground.  Only  one  pilot  problem  has  been  run  and  the  coiqputed  results  are  ex¬ 
tremely  Interesting  although  Inadequate  for  use  as  design  Inputs.  Certain  nec¬ 
essary  assumptions  were  found  to  be  unrealistically  chosen;  this  caused  some 
Inconsistencies.  Such  occurrences  are,  of  course,  the  rule  for  first  results 
from  highly  complex  analyses,  and  do  not  affect  the  validity  of  the  procedure. 

The  preliminary  results,  which  Include  the  effect  of  air-blast  loading  on  the 
surface,  show  large  attenuation  of  stress  euod  energy  in  the  plastic  zone.  Mo 
resvilts  were  obtained  near  the  surface,  however,  because  the  shallowest  ray 
which  was  used  was  at  an  angle  of  l8  degrees  from  the  horizontal.  The  solution 
Is  capable  of  considering  shallower  rays  (even  surface  rays)  and  such  rays  should 
be  Included  In  any  extensions  of  the  computatlcn  s .  The  most  Important  short¬ 
coming  of  the  computer  program,  from  which  these  results  have  been  obtained,  Is 
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its  low  speed-  Another  lack  is  that  of  a  more  sophisltlcated  yield  condition 
than  that  used  thus  far.  After  these  probleos  are  corrected,  paraaetric  studies 
should  be  undertaken  to  show  how  sensitive  results  are  to  ground  aaterlal  prop¬ 
erties.  When  this  has  been  established,  typical  earth  naterials  properties 
should  be  used  in  confutations  \diich  would  be  suitable  for  design  inputs  for 
structures  which  might  be  located  in  the  plastic  medium- 

In  addition,  based  on  Sewmaxk's  conclusion  that  Interface  discrepancies  do  not 
seriously  affect  the  pbenomsna  transmitted  across  the  Interface,  it  will  be  pos¬ 
sible  to  cotrislne  the  hydrodynamic,  plastic,  and  elastic  analyses  in  one  solution, 
and,  further,  to  include  alr-lndueed  effects-  nils  conbinatlon  will  constitute 
the  iirst  cosflete  and  unified  description  of  ground  effects  in  the  close-in 
region  and  beyond  Proposed  means  of  accoopllshing  this  are  outlined  in  the 
following  section . 

B.  PROPOSED  METHODS  FOR  CCMBIXIBQ  EFFECTS. 

The  Brooks  analysis  has  been  shown  to  be  capable  of  handling  the  ground  stresses 
and  motions  associated  with  the  zone  of  hydrodynamic  behavior.  The  precise  boun¬ 
dary  at  \diich  behavior  changes  from  hydrodynamic  to  plastic  is  unknown  and, 
indeed,  probably  does  not  exist  as  a  definable  interface.  Therefore,  a  suitable 
radius  can  be  chosen,  as  was  done  for  the  Brooks  pilot  problem,  and  Inputs  to  the 
plastic  analysis  obtained  from  the  hydrodynamic  stresses  euid  motions.  In  the 
hydrodynamic  region  air  blast  effects  are  trivial.,  although  the  hydrodynamic 
analysis  does  include  such  effects. 
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nw  plastic  analysis^  Including  air  blast  effects  (v/hich  Incrcucc  In  in^xir- 
tance  with  Increasing  distance  from  the  burst)  gives  results  which  can  be 
used  to  extend  the  description  of  ‘Uie  phenowna  into  the  elastic  regime. 

This  can  be  accomplished  in  three  possible  ways. 

The  first  method  is  to  continue  the  elasto-plastic  analysis  according  to 
the  Brooks  equation  into  the  elastic  zone.  Ihls  would  not  result  In  the  most 
detailed  estimates  of  stresses  and  strains  because  the  procedure  Is  hindered 
by  the  non-linear  features  of  the  plastic  analysis  which  are  not  required  for 
the  elastic  solution. 

The  second  method  is  to  use  the  Newmeirk  cuialysls  in  idiich  the  stress  inputs 
Into  the  elastic  zone  are  obtained  from  the  Brooks  plastic  analysis.  This 
procedure  presents  no  intractable  difficulty  and  appears  rather  attractive. 

The  Mewmark  euialysls  nusti  of  course^  be  extended  to  points  nearer  the  sur¬ 
face  at  distances  of  Interest  and  air  blast  effects  must  be  added  to  the  des¬ 
cription  ■ 

The  last  method  is  to  a  half-space  solution  similar  to  that  of  freeman 
Gilbert  for  the  Rayleigh  wave  problem.  The  procedure  would  consist  of 
applying  to  the  surface  of  the  elastic  solid  a  point  source  Input  \diich  is 
derive^  from  the  energy  attenuation  factor  computed  in  the  plastic  analysis. 
This  procedure  Is  feasible  because  no  energy  attenuation  takes  place  In  an 
elastic  body.  The  characteristics  of  the  surface  Input  can  be  varied  until 
a  match  is  obtained  between  the  stresses  and  motion  (at  the  radius  corres¬ 
ponding  to  the  plastic  zone  limit)  in  the  elastic  half- space  and  those  obtained 
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from  the  confined  hydrodynamic  and  plastic  analyses.  It  will  be  necessary 
to  choose  the  source  such  that  air  blast  effects  are  the  as  those  from  a 
surface  burst  or,  alternatively,  to  use  superposition  to  coobine  a  solution 
for  an  air  burst  (such  as  the  Qilbert  Bayleigh  vave  pulse)  with  the  surface 
source.  Since  the  solutions  are  llnear,no  difficulties  arise.  This  last  ap¬ 
proach  has  certain  advantages  in  dealing  with  locations  near  the  surface  but 
is  not  as  versatile  as  Bewaark's  numerical  technique. 
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